Photon storage in A-type optically dense atomic media. II. Free-space model 
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In a recent paper [Gorshkov et al, Phys. Rev. Lett. 98, 123601 (2007)], we presented a universal 
physical picture for describing a wide range of techniques for storage and retrieval of photon wave 
packets in A-type atomic media in free space, including the adiabatic reduction of the photon 
group velocity, pulse-propagation control via off-resonant Raman techniques, and photon-echo based 
techniques. This universal picture produced an optimal control strategy for photon storage and 
retrieval applicable to all approaches and yielded identical maximum efficiencies for all of them. In 
the present paper, we present the full details of this analysis as well some of its extensions, including 
the discussion of the effects of non-degeneracy of the two lower levels of the A system. The analysis 
in the present paper is based on the intuition obtained from the study of photon storage in the 
cavity model in the preceding paper [Gorshkov et al., Phys. Rev. A 76, 033804 (2007)]. 
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I. INTRODUCTION 

High fidelity storage of a traveling pulse of light into 
an atomic memory and the subsequent retrieval of the 
state back onto a light pulse are currently being pur- 
sued by a number of laboratories around the world. In 
a recent paper we presented a universal picture for 
describing, optimizing, and showing a certain degree of 
equivalence between a wide range of techniques for pho- 
ton storage and retrieval in A-type atomic media, includ- 
ing the approaches based on electromagnetically induced 
transparency (EIT), off-resonant Raman processes, and 
photon echo. In the present paper, as well as in the pre- 
ceding paper (which wc will refer to henceforth as 
paper I) and the paper that follows (which we will 
refer to henceforth as paper III), we present all the de- 
tails behind this universal picture and the optimal control 
shaping that it implies, as well as consider several exten- 
sions of this analysis beyond the results of Ref. In 
particular, in paper I, we present the full details of the op- 
timization in the slightly simpler model where the atoms 
are placed inside a cavity. Using the intuition gained 
from the cavity model discussion, we show in the present 
paper all the details behind the analysis of the free-space 
model given in Ref. We also discuss several exten- 
sions of the analysis of Ref. [l[ , such as the inclusion of 
the decay of coherence between the two lower levels of the 
A system and the effects of nondegeneracy of these two 
levels. Finally, in paper III, we generalize our treatment 
to include the effects of inhomogeneous broadening. 

For a complete introduction to photon storage in A- 
type atomic media, as it applies to paper I and to the 
present paper, as well as for the full list of references, we 
refer the reader to paper I. In the present Introduction, 
we only list the two main results of the present paper. 
The first important result is the abovementioned proof 
of a certain degree of equivalence between a variety of 
different photon storage protocols. In particular, this 



result means that provided there is a sufficient degree 
of control over the shape of the incoming photon wave 
packet and/or over the power and shape of the classical 
control pulses, all the protocols considered have the same 
maximum achievable efficiency that depends only on the 
optical depth d of the medium. The second important re- 
sult is a novel time-reversal-based iterative algorithm for 
optimizing quantum state mappings, a procedure that 
we expect to be applicable beyond the field of photon 
storage. One of the key features of this optimization al- 
gorithm is that it can not only be used as a mathematical 
tool but also as an experimental technique. In fact, fol- 
lowing our theoretical proposal, an experimental demon- 
stration of this technique has already been carried out 
[U • Both the experimental results [1] and the theoretical 
results of the present paper indicate that the suggested 
optimization with respect to the shape of the incoming 
photon wave packet and/or the control pulse shape and 
power will be important for increasing the photon storage 
efficiencies in current experiments. 

Although the slightly simpler cavity model discussed 
in paper I is similar enough to the free-space model to 
provide good intuition for it, the two physical systems 
have their own advantages and disadvantages, which we 
will discuss in the present paper. One advantage of the 
free-space model is the fact that it is easier to set up ex- 
perimentally, which is one of the reasons we study this 
model in the present paper. Turning to the physics of 
the two models, the main differences come from the fact 
that in the cavity model the only spin wave mode acces- 
sible is the one that has the excitation distributed uni- 
formly over all the atoms. In contrast, in the free-space 
model, incoming light can couple to any mode specified 
by a smooth excitation with position-dependent ampli- 
tude and phase. As a consequence of this, the free-space 
model allows for high efficiency storage of a wider range 
of input light modes than the cavity model. In particu- 
lar, we showed in paper I that in the cavity model high 
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efficiency photon-echo-based storage (which we refer to 
as fast storage) is possible for a single input mode of 
duration ~ 1/(7(7), where 7 is the optical polarization 
decay and C is the cavity cooperativity parameter. In 
contrast, we show in the present paper that high effi- 
ciency fast storage in a free-space atomic ensemble with 
optical depth d is possible for any input light mode of 
duration T provided T7 <C 1 and Tdj 3> 1. However, 
the cavity model also has some advantages over the free- 
space model. In particular, the error during optimal light 
storage and retrieval for a given atomic ensemble scales 
as the inverse of the optical depth, as we have shown for 
the cavity model in paper II and for the free-space model 
in the present paper. The optimal efficiency is therefore 
higher when the ensemble is enclosed in a cavity, which 
effectively enhances the optical depth by the cavity fi- 
nesse to form the cooperativity parameter C. Moreover, 
if one is forced to retrieve from a spatially uniform spin 
wave mode (e.g., if the spin wave is generated via spon- 
taneous Raman scattering Q), the error during retrieval 
will decrease faster with optical depth in the cavity model 
(~ 1/C 1 ) than in the free-space model (~ 1/Vd), 

The remainder of the present paper is organized as fol- 
lows. In Sec. HH the model is introduced. In Sees. IIII| 
IIV1 and [V] we prove that during retrieval there exists a 
fixed branching ratio between the desired light emission 
rate and undesired polarization decay rate, and use this 
in combination with time reversal to derive the optimal 
strategy for storage and retrieval without fully solving the 
equations. In Sees. IVIIand IVII1 the equations are solved 
analytically in the adiabatic and fast limits, respectively, 
and more specific statements about the optimal control 
strategy are made. In Sec. IVIII1 the effect of nondegener- 
acy of the two metastable states is discussed. In Sec. IIX1 
we summarize the discussion of the free-space model. Fi- 
nally, in the Appendixes, we present some details omitted 
in the main text. 



II. MODEL 

We refer the reader to AppendixlAlfor the details of the 
model and for the derivation of the equations of motion. 
In this section, we only briefly summarize the model and 
state the equations of motion without derivation. 

We consider a free-space medium of length L and cross- 
section area A containing N = J Q dzn(z) atoms, where 
n(z) is the number of atoms per unit length. We assume 
that within the interaction volume the concentration of 
atoms is uniform in the transverse direction. The atoms 
have the same A-type level configuration as in the cavity 
case discussed in paper I and shown in Fig. 1 of paper I. 
They are coupled to a quantum field and a copropagating 
classical field. We assume that quantum electromagnetic 
field modes with a single transverse profile are excited. 
We also assume that both the quantum and the classi- 
cal field are narrowband fields centered at uj\ = cu eg — A 
and ui2 = uj es — A, respectively (where u> eg and oj es are 



atomic transition frequencies). The quantum field is de- 
scribed by a slowly varying operator £(z,t), while the 
classical field is described by the Rabi frequency enve- 
lope £l(z, t) = Q,(t — z/c). 

We neglect reabsorption of spontaneously emitted pho- 
tons. This is a good approximation since we are inter- 
ested in the storage of single- or few-photon pulses, in 
which case there will be at most a few spontaneously 
emitted photons. Although for an optically thick medium 
they can be reabsorbed and reemitted jy, 0] , the proba- 
bility of spontaneously emitting into the mode £ is given 
by the corresponding far-field solid angle ~ A 2 /A ~ d/N, 
where A is the cross section area of both the quantum 
field mode and the atomic medium (see Appendix [X] 
for a discussion of why this choice is not important), 
A = 2-kc/uj! is the wavelength of the quantum field, and 
d ~ A 2 N/A is the resonant optical depth of the ensem- 
ble. In most experiments, this probability is very small. 
Moreover, we will show that for the optimized storage 
process, the fraction of the incoming photons lost to 
spontaneous emission will decrease with increasing op- 
tical depth. In practice, however, reabsorption of spon- 
taneously emitted photons can cause problems [f| dur- 
ing the optical pumping process, which is used to initial- 
ize the sample, and this may require modification of the 
present model. 

We treat the problem in a one-dimensional approx- 
imation. This is a good approximation provided that 
the control beam is much wider than the single mode 
of the quantum field defined by the optics, as, for ex- 
ample, in the experiment of Ref. [f|. In this case, the 
transverse profile of the control field can be considered 
constant; and, in the paraxial approximation, the equa- 
tions reduce to one-dimensional equations for a single 
Hermite-Gaussian quantum field mode 0, [Toj . 

We define the polarization operator P(z,t) = 
\Na ge (z,t) and the spin-wave operator S(z,t) = 
Va" ijgs (z, t) (where a^ v (z, t) are slowly varying position- 
dependent collective atomic operators defined in Ap- 
pendix [A|. In the dipole and rotating- wave approxima- 
tions, to first order in £ , and assuming that at all times 
almost all atoms are in the ground state, the Heisenberg 
equations of motion read 

(d t + cd z )£ = igv f NPn(z)L/N, (1) 
dtP = -(i + iA)P+igVN£+inS+^Fp, (2) 
d t S= - ls S + ifl*P+ y^sFs, (3) 

where we introduced the spin- wave decay rate j s , the po- 
larization decay rate 7, and the corresponding Langevin 
noise operators Fp(z,t) and Fs(z,t). As in the cavity 
case, collective enhancement [1J| results in the increase 
of the atom-field coupling constant g (assumed to be real 
for simplicity) by a factor of VN up to gVN. 

As we show in Appendix [A] and explain in detail in 
paper I, under reasonable experimental conditions, the 
normally ordered noise correlations of Fp and F$ are 
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zero, i.e., the incoming noise is vacuum and the trans- 
formation is passive. As we show in Sec. II of paper I, 
this implies that efficiency is the only number required 
to completely characterize the mapping. 

As in the cavity discussion of paper I, we suppose that 
initially all atoms are in the ground state, i.e., no atomic 
excitations are present. We also assume that there is only 
one nonempty mode of the incoming quantum field and 
that it has an envelope shape ho(t) nonzero on [0, T] . The 
term "photon storage and retrieval" refers to mapping 
this mode onto some mode of 5 and, starting at a later 
time T T > T, retrieving it onto an outgoing field mode. 
Then precisely as in the cavity case in paper I, for the 
purposes of finding the storage efficiency, which is given 
by the ratio of the numbered of stored excitations to the 
number of incoming photons 



J L dz^{S\ Zl T)S{z,T)) 



(4) 



we can ignore Fp and F$ in Eqs. CO)-® and treat these 
equations as complex number equations with the inter- 
pretation that the complex number fields describe the 
shapes of quantum modes. In fact, although the result- 
ing equations describe our case of quantized light coupled 
to the \g) — |e) transition, they will also precisely be the 
equations describing the propagation of a classical probe 
pulse. To sec this, one can simply take the expectation 
values of Eqs. CD)-© an d use the fact that classical probe 
pulses are described by coherent states. 

Going into a comoving frame t' = t — z/c, introducing 
dimensionless time t = jt' and dimensionless rescaled 
coordinate z — J Q dz'n(z')/N, absorbing a factor of 

^Jc/(L^) into the definition of £ , we reduce Eqs. ©-© 
to 



d~ z £ = iVdP, 

d~ t P = -(l+iA)P + iVd£ 
fyS = ih*(i)P, 



(5) 
(6) 
(7) 



where we have identified the optical depth d = 
g 2 NL/( r yc) and where A = A/7 and = fi/7. We 
confirm in Appendix |A1 that from the definition it follows 
that d is independent of the size of the beam and, for a 
given transition, only depends on the density and length 
of the ensemble. Moreover, the definition of d that we 
use here can be related to the intensity attenuation of a 
resonant probe in our three level system with the control 
off, in which case the equations give an attenuation of 
exp(— 2c?). In Eqs. I©-© and in the rest of this paper 
(except for Sec. I VIE]) , we neglect the decay 7 S of the spin 
wave. However, precisely as in the cavity case, nonzero 
7 S simply introduces an exponential decay without mak- 
ing the solution or the optimization harder, as we will 
discuss in Sec. I VI El We also note that Eqs. CO)-© are 
the same as the equations of Ref. [12] for copropagating 
fields, generalized to nonzero A and 7 S , and taken to first 
order in £ . 



During storage, shown (in original variables) in 
Fig. [IJa), the initial and boundary conditions are (in 
rescaled variables) £(z = 0, t) = £- m {t), P(z,i = 0) = 0, 
and S(z, t = Q) = 0, where £i n (t) is nonzero for i S [0, T] 
(where T — T7) and, being a shape of a mode, is nor- 
malized according to / Q T di\£- m {i)\ 2 = 1. S{z,T) gives the 
shape of the spin-wave mode, into which we store, and 
the storage efficiency is given by 



dz\S{z,f){ 



(8) 



Loss during storage comes from the decay 7 as well as 
from the "leak" £(z — l,t) shown in Fig. [TJa). Then 
at a later time T r > T (where T r — T r -y), we want to 
retrieve the excitation back onto a photonic mode cither 
in the forward direction, as shown in Fig. [ljb), or in 
the backward direction [13| (i.e., with the retrieval con- 
trol pulse incident from the right) as shown in Fig.QJc). 
Instead of turning our Eqs. ©-© around to describe 
backward retrieval, we invert, for backward retrieval, the 
spin wave according to S(z,T r ) = 5(1 — z,T), whereas 
we keep S(z, T T ) = S(z, T) for forward retrieval. Because 
of the z-dependent phases in Eq. (|A10[) , this prescrip- 
tion for backward retrieval is strictly valid only for zero 
splitting between the two metastable states (oj sg = 0). 
In Sec. IVIII1 we will discuss the effect of nonzero oj sg . 
The remaining initial and boundary conditions during 
retrieval are £(z = 0,t) = and P(z,T r ) = 0. If we 
renormalize the spin wave before doing the retrieval, then 
the retrieval efficiency will be given by 



dt\£(l,i)\ 7 



(9) 



If we do not renormalize the spin wave before doing the 
retrieval, this formula will give the total efficiency of stor- 
age followed by retrieval r^ot = ?7s?7r- 

To solve Eqs. ©-©, it is convenient to Laplace trans- 
form them in space according to z — > u, so that Eqs. do I 
become 



£ 



■ Vd £ in 
1 P H 

u u 



(10) 



dfP = -{l + - + iA)P + tn{t)S + i—£ iD . (11) 
u u 

As in the cavity case in paper I, it is also convenient to 
reduce Eqs. ©-© to a single equation 



5-^5 



-{i + -+iA)s+\n\ 2 s = -n*—£ in , (12) 

u u 



where the overdot stands for the t derivative. As in the 
cavity case, this second-order differential equation can- 
not, in general, be fully solved analytically. Similar to 
paper I, we can, however, derive several important re- 
sults regarding the optimal control strategy for storage 
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Sin(0 




FIG. 1: (Color online) (a) Storage, (b) forward retrieval, and 
(c) backward retrieval setup. The smooth solid curve is the 
generic control field shape (Q) for adiabatic storage or re- 
trieval; the dotted square pulse indicates a 7r-pulse control 
field for fast storage or retrieval. The dashed line indicates the 
quantum field £ and the spin-wave mode S. During storage, 
£(L,t) is the "leak," whereas it is the retrieved field during 
retrieval. 



and retrieval without making any more approximations. 
We discuss these results in Sec. IIII) where we optimize 
retrieval, and in Sees. IIVI and W\ where we introduce the 
important time reversal ideas, which allow us to deduce 
the optimal storage from the optimal retrieval. 



III. OPTIMAL RETRIEVAL 

Although Eq. (JT2J) cannot, in general, be fully solved 
analytically, we still make in this and in the following two 
sections several important statements regarding the op- 
timal strategy for maximizing the storage efficiency, the 
retrieval efficiency, and the combined (storage followed 
by retrieval) efficiency without making any more approx- 
imations, ft is convenient to first consider retrieval, and 
we do so in this section. 

Although we cannot, in general, analytically solve for 
the output field £ ut(t), we will show now that, as in the 
cavity case in paper f, the retrieval efficiency is indepen- 
dent of the control shape and the detuning provided no 



excitations are left in the atoms. Moreover, we will show 
that the retrieval efficiency is given by a simple formula 
that depends only on the optical depth and the spin- wave 
mode. From Eqs. ([7]) and (fTTj) . it follows that 



% (P(u, t) [P(u'*,t)j * + S(u, t) [S(u'*, t)] ' 

= -{2 + d/u + d/u')P(u,t) [P(u'*,t)]* . (13) 

Using Eqs. (fTU|) and (fT3"|) and assuming P(u,oo) = 
S(u, oo) = (i.e., that no excitations are left in the atoms 
at t = oo), the retrieval efficiency is 



Vr = £- 1 — / dtP(u, t) [P{u'* , t) 



c- 1 



—p-S[u, T r ) S{u'\T t ) 



2uu' + d(u + u') 

dz / dz'S(l - z,T,:)S*(l - z',T T )kr(z,z'),{U) 
o Jo 

where L~ l means that two inverse Laplace transforms 
(u — > z and u' — > z') are taken and are both evaluated at 
z = z! = 1 and where the kernel k r is defined as 



k r (z, z')=L 



2uu' + d(u + u') 



i(l — z)—u'(l — z') 



-e-d^IoidVUP), 



(15) 



where /„ is the nth-order modified Bessel function of the 
first kind. Similarly to the cavity case in paper 1, we see 
that the efficiency is independent of A and f2, which re- 
flects that in Eq. (fTTj) (or, equivalently, on the right-hand 
side of Eq. (fT3|) ) there is a fixed branching ratio between 
the decay rates of P. For a given u the rates are (in the 
original units) 7 and ^d/u into the undesired modes and 
the desired mode £ ut, respectively, independent of A 
and Q. In fact, a stronger result than the independence 
of retrieval efficiency from A and Q can be obtained: as 
we show in Appendix [B] the distribution of spontaneous 
emission loss as a function of position is also independent 
of the control and detuning. 

In contrast to the cavity case in paper I where there 
was only one spin-wave mode available, in the free-space 
case, the retrieval efficiency in Eq. (fT4"|) is different for 
different spin-wave modes. We can, thus, at each d, opti- 
mize retrieval by finding the optimal retrieval spin wave 
Sd{z) (we suppress here the argument T r ). The expres- 
sion for the efficiency in the last line of Eq. (|14[) is an 
expectation value of a real symmetric (and hence Hermi- 
tian) operator k T (z, z') in the state S'(l — z). It is there- 
fore maximized when S(l — z) is the eigenvector with the 
largest eigenvalue of the following eigenvalue problem: 



VvS(l - z) 



dz'k r {z,z')S(l - z ). 



(16) 



Since eigenvectors of real symmetric matrices can be cho- 
sen real, the resulting optimal spin wave Sd{z) can be 
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chosen real. To find it, we start with a trial S(z) and 
iterate the integral in Eq. (TTo| several times until con- 
vergence In Fig. [21 we plot the resulting optimal 
spin wave Sd(z) for d = 0.1, 1, 10, 100, as well as its lim- 
iting shape (5oo(^) = V^z) as d — ► oo. At d — > 0, the 
optimal mode approaches a flat mode. These shapes can 
be understood by noting that retrieval is essentially an 
interference effect resembling superradiance, where the 
emission from all atoms contributes coherently in the for- 
ward direction. To get the maximum constructive inter- 
ference, it is desirable that all atoms carry equal weight 
and phase in the spin wave. In particular, at low opti- 
cal depth, this favors the flat spin wave. On the other 
hand, it is also desirable not to have a sudden change 
in the spin wave (except near the output end of the 
ensemble). The argument above essentially shows that 
excitations can decay through two different paths: by 
spontaneous emission in all directions or by collective 
emission into the forward direction. In Eq. ©, these 
two paths are represented by the — P and i\fd£ terms, 
respectively. The latter gives rise to a decay because 
Eq. ([5]) can be integrated to give a term proportional to 
P: £ — i J dzydP. To obtain the largest decay in the 
forward direction all atoms should ideally be in phase so 
that the phase of P(z) is the same at all z. This construc- 
tive interference, however, is not homogeneous but builds 
up through the sample. At z = 0, we have £ = 0, and 
the spontaneous emission is, thus, the only decay chan- 
nel, i.e., d\P(z = 0)\ 2 /dt = -2\P(z = 0)| 2 . To achieve 
the largest retrieval efficiency, we should therefore put a 
limited amount of the excitation near z = and only 
have a slow build up of the spin wave from z = to 
2 = 1. The optimal spin- wave modes in Fig. [5] represent 
the optimal version of this slow build up. We will also 
reinterpret these optimal modes from a different perspec- 
tive in Sec. IVTAl using the EIT window concept. 

From the qualitative argument given here, one can es- 
timate the dependence of the optimal retrieval efficiency 
on the optical depth: we consider the emission into a for- 
ward mode of cross sectional area A. In the far field, this 
corresponds to a field occupying a solid angle of \ 2 /A, 
where A = (2irc)/u>i is the wavelength of the carrier. A 
single atom will, thus, decay into this mode with a proba- 
bility ~ A 2 1 A. With N atoms contributing coherently in 
the forward direction, the emission rate is increased by a 
factor of AT to 7f ~ 7AA 2 /A The retrieval efficiency can 
then be found from the rate of desired (tj) and undesired 
(7) decays as r) = 7f/(7f + 7) ~ 1 - A/ {NX 2 ). By not- 
ing that A 2 is the cross section for resonant absorption 
of a two-level atom, we recognize AA 2 /^4 as the optical 
depth d (up to a factor of order 1). The efficiency is then 
77 ~ 1 — l/rf, which is in qualitative agreement with the 
results of the full optimization which gives 1 — 77 ss 2.9 jd. 
A more detailed discussion of the dependence of retrieval 
efficiency on the shape of the spin wave is postponed un- 
til Sec. lVI Al where Eq. (TH)l is rederived in the adiabatic 
limit and discussed in the context of the EIT window. 



S d (z) 




FIG. 2: Optimal modes §d(£) to retrieve from (in the forward 
direction) at indicated values of d. The flipped versions of 
these modes Sd(l — z) are the optimal modes for backward 
retrieval and are also the optimal (normalized) spin waves 
S(z , T) I -^/77f lax for adiabatic and fast storage if it is optimized 
for storage alone or for storage followed by backward retrieval 
(77™ ax is the maximum storage efficiency). 



IV. OPTIMAL STORAGE FROM TIME 
REVERSAL: GENERAL PROOF 

As already mentioned in paper I, the concept of time 
reversal allows us to deduce the optimal control strategy 
for storage from retrieval. In this section, we prove this 
result both for the free-space case and for the cavity case 
(the cavity case differs only in that there is just one spin- 
wave mode involved). In the next section, we generalize 
these ideas and show that time reversal can be generally 
used to optimize state mappings. 

Despite the fact that our system contains nonreversible 
decay 7, time reversal is still an important and meaning- 
ful concept: to make time reversal applicable in this situ- 
ation, we expand our system so that we not only consider 
the electric field and the spin wave, but also include all 
the reservoir modes, into which the excitations may de- 
cay. As discussed in Sec. Inland Appendix [Al the initial 
state of the reservoir modes is vacuum. When consider- 
ing "all the modes in the universe" [l5| , we have a closed 
system described by (possibly infinitely many) bosonic 
creation operators {Oj} with commutation relations 



o l ,o] 



= Si 



(17) 



The evolution we consider here can be seen as a gener- 
alized beam-splitter transformation, and can equivalently 
be specified as a Heisenberg picture map between the 
annihilation operators (\out = Ylj Uij[T,0; Ct(t)]Oj t i a 

or as a Schrodinger picture map U[T, 0; = 
J2ij Uij[T,0;Q(t)]\i){j\ in the Hilbert space U with 
an orthonormal basis of single excitation states \i) — 
0\ I vacuum). To stress that the mapping depends on 
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the classical control field 0(t), we here include the ar- 
gument Q(t) in the evolution operator U[t 2 , n; 
which takes the state from time n to t%. The op- 
erator U[T, 0; Q(t)] must be unitary U*[T, 0; Q(t)] = 
U^iT^iQit)] = U[0,T;Q(t)]. For simplicity of nota- 
tion, we will here use the Schrodinger picture. 

Let us define two subspaces of TL: subspace A of "ini- 
tial" states and subspace B of "final states." B , the 
orthogonal complement of B, can be thought of as the 
subspace of "decay" modes (that is, the reservoir and 
other states, possibly including the "initial" states, to 
which we do not want the initial states to be mapped). 
In this section, we will use U as the retrieval map, in 
which case A and B are spin-wave modes and output 
photon modes, respectively, while B ± includes A, empty 
input field modes, and the reservoir modes, to which the 
excitations can decay by spontaneous emission. 

In the cavity derivation in paper I, we solved in the 
adiabatic limit for the control pulse shape f2 r (i), which 
retrieves the atomic excitation into a specific mode e(t). 
We then derived the pulse shape f2 s (i), which optimally 
stores an incoming mode £ m (t) , and noted that if the in- 
coming mode is the time-reverse of the mode, onto which 
we retrieved, i.e., £* n (T—t) — e(t), then the optimal stor- 
age control is the time-reverse of the retrieval control, i.e., 
f2*(T — t) — Sl T (t). Furthermore, in this case the storage 
and retrieval efficiencies were identical. As we now show, 
this is not a coincidence, but a very general result. 

For the free-space case, we define the "overlap effi- 
ciency" for storing into any given mode S(z) as the ex- 
pectation value of the number of excitations in this mode 
S(z). Since the actual mode (call it S'(z)), onto which 
the excitation is stored, may contain components orthog- 
onal to S(z), the overlap efficiency for storing into S(z) 
is in general less than the (actual) storage efficiency, and 
is equal to it only if S(z) = S'(z). 

We will now prove that storing the time reverse of the 
output of backward retrieval from S* (z) with the time re- 
verse of the retrieval control field gives the overlap storage 
efficiency into S(z) equal to the retrieval efficiency. To 
begin the proof, we note that the probability to convert 
under U an initial excitation from a state \a) in A into a 
state | b) in B is just 

r 1 =\{b\U[T 1 Q-Mma)\ 2 = MU-\T^Mt)]\b)\\ (18) 

where, in the last expression, we have used the unitar- 
ity of U[T,0;Cl(t)]. We now assume that U[T,Q;Q(t)] 
describes retrieval and that \a) stands for S*(z), while 
| b) stands for the output field mode £, onto which S* 
is retrieved under f2(£). Then r\ is just the retrieval effi- 
ciency from S* . The last expression then shows that if 
we could physically realize the operation J7 -1 [T, 0; f2(t)], 
then it would give the overlap efficiency for storage of 
£ into S* equal to the retrieval efficiency ij. The chal- 
lenge is, therefore, to physically invert the evolution and 
realize [7 — 1 [T, 0; Q(t)]. As we now show, time reversal 
symmetry allows us to perform this inverse evolution in 



some cases. We refer the reader to Appendix [C] for a 
careful definition of the time reversal operator T and for 
the proof of the following equality: 

U-^T, 0; ft(t)} = TU[T, 0; fT(T - t)]t, (19) 

where it is implicit that the carrier wave vector of the 
time-reversed control pulse 57* (T — f) propagates in the 
direction opposite to the carrier of f2(i). Physically, 
Eq. (|19p means that we can realize the inverse evolu- 
tion by time-reversing the initial state, evolving it using 
a time-reversed control pulse, and finally time-reversing 
the final state. Then, using Eqs. (fT5| and (fT!)|) , the re- 
trieval efficiency may be rewritten as 

V = \(a\fU[T,0;^(T-t)]f\b)\ 2 . (20) 

This means that if we can retrieve the spin wave S* back- 
wards onto £{t) using we can use fi*(T — t) to store 
£m{t) — £*{T — t) onto S with the overlap storage effi- 
ciency equal to the retrieval efficiency 

We will now prove that this time-reversed storage is 
also the optimal solution, i.e., that an overlap efficiency 
for storage into <S* greater than rj r is not possible. To be- 
gin the proof, let us suppose, on the contrary, that we can 
store &m(t) into S with an overlap efficiency rj s > r] T . Ap- 
plying now the time reversal argument to storage, we find 
that backward retrieval from S* with the time-reversed 
storage control will have efficiency greater than r\ x . How- 
ever, from Eq. (|14[) . we know that the retrieval efficiency 
is independent of the control field and is invariant under 
the complex conjugation of the spin wave, so we have 
reached a contradiction. Therefore, the maximum over- 
lap efficiency for storage into a given mode S is equal to 
the backward retrieval efficiency from S* (and S) and can 
be achieved by time-reversing backward retrieval from 
S*. 

Finally, the strategy for storing £- m {t) with the max- 
imum storage efficiency (rather than maximum overlap 
efficiency into a given mode, as in the previous para- 
graph) follows immediately from the arguments above: 
provided we can retrieve the (real) optimal backward- 
retrieval mode §d(L — z) backwards into £* n (T — t), the 
optimal storage of £- m if) will be the time reverse of this 
retrieval and will have the same efficiency as the optimal 
retrieval efficiency at this d, i.e., the retrieval efficiency 
from Sd- 

While the above argument is very general, it is impor- 
tant to realize its key limitation. The argument shows 
that it is possible to optimally store a field £i a (t) provided 
we can optimally retrieve onto £* n (T — t) (i.e., backward- 
retrieve Sd(L — z) into £* n (T — t)). It may, however, not 
be possible to optimally retrieve onto £^{T — t) because 
it may, for example, be varying too fast. For this reason, 
we shall explore in Sees. I VI Al and FyTA onto which fields 
it is possible to retrieve a given spin wave. Before we 
do this, however, we will show in the next section that 
time reversal does not only allow one to derive the opti- 
mal storage strategy from the optimal retrieval strategy, 
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as we did in this section, but also allows one to find the 
optimal spin wave for retrieval. 



V. TIME REVERSAL AS A TOOL FOR 
OPTIMIZING QUANTUM STATE MAPPINGS 

We will now show that time reversal can be used as 
a general tool for optimizing state mappings. Moreover, 
we will show that for the photon storage problem consid- 
ered in this paper, in addition to being a very convenient 
mathematical tool, the optimization procedure based on 
time reversal may also be realized experimentally in a 
straightforward way. 

In Sec. IIIII we found Sd(z), the optimal spin wave to 
retrieve from, by starting with a trial spin wave S%(z) 
and iterating Eq. (|16|) until convergence. While we just 
used this as a mathematical tool for solving an equation, 
the iteration procedure actually has a physical interpre- 
tation. Suppose that we choose a certain classical control 
fi(t) and retrieve the spin wave S(z) forward onto £ (t) 
and then time reverse the control to store £*{T — t) back- 
wards. By the argument in the last section, this will, in 
general, store into a different mode S'(z) with a higher 
efficiency (since the actual storage efficiency is, in gen- 
eral, greater than the overlap storage efficiency into a 
given mode). In this way, we can iterate this procedure 
to compute spin waves with higher and higher forward 
retrieval efficiencies [HI Il7j . In fact, forward retrieval 
followed by time-reversed backward storage can be ex- 
pressed as 

S 2 {l-z)= [ d?Jk I (S,S')S* 1 (l- z'), (21) 
Jo 

which for real S is identical to the iteration of Eq. (JTBJ). 
We note that the reason why backward storage had to 
be brought up here (in contrast to the rest of the paper, 
where storage is always considered in the forward direc- 
tion) is because Eq. ([T5)| . which Eq. (|2"Tj) is equivalent to 
for real S, discusses forward retrieval, whose time-reverse 
is backward storage. 

Since the iterations used to maximize the efficiency in 
Eq. (TT(3|) are identical to Eq. (f2"Tj) . the physical interpre- 
tation of the iterations in Eq. (fl6|) is that we retrieve 
the spin wave and store its time-reverse with the time- 
reversed control field (i.e., implement the inverse C/ _1 of 
the retrieval map U using Eq. §T§\i ). We will explain 
below that this procedure of retrieval followed by time- 
reversed storage can be described mathematically by the 
operator /';/' • where Pa and Pb are the pro- 

jection operators on the subspaces A and B of spin wave 
modes and output photon modes, respectively, and where 
Af provides renormalization to a unit vector. It is, in fact, 
generally true that in order to find the unit vector \a) 
in a given subspace A of "initial" states that maximizes 
the efficiency i] — |Pg£/|a)| 2 of a given unitary map U 
(where B is a given subspace of "final" states), one can 



start with any unit vector \a) G A and repeatedly apply 
NPaU~ 1 PbU to it. We prove in Appendix |D] that this 
procedure converges to the desired optimal input mode 
l^max) yielding the maximum efficiency ?7 m ax (provided 
(a|a max ) 7^ 0) and that MpBU\a ma x) also optimizes U^ 1 
as a map from B to A. 

We have just discussed the iterative optimization pro- 
cedure as a purely mathematical tool for computing the 
optimal initial mode in the subspace A of initial states. 
However, in the previous section and in Appendix [CJ we 
showed that for our system one can implement the inverse 
evolution U^ 1 experimentally by first time reversing the 
state of the system, then applying ordinary evolution but 
with a time-reversed control, and then time reversing the 
state of the system again. Thus, in addition to being a 
convenient mathematical tool, the time reversal based 
optimization technique suggested in this paper is an ex- 
perimentally realizable procedure. In fact, following this 
work, this procedure has recently been applied experi- 
mentally to the optimization of the storage and retrieval 
of light [|. 

As an example, let us discuss how this experimental 
implementation applies to the optimization of backward 
retrieval. The full Hilbert space is spanned by subspace A 
of spin-wave modes, subspace B of output field modes, as 
well as a subspace containing (empty) input and reservoir 
field modes. The goal is to optimize the retrieval time 
evolution map E/[0(t)] for some fixed detuning A and 
fixed backward-propagating Rabi frequency pulse fl(t) 
(sufficiently powerful for complete retrieval) with respect 
to the initial spin wave \a) € A. From Sec. IIV1 it follows 
that the iteration AfP A U- 1 [n(t)]pBU[^l(t)]\a) required 
for the optimization can be experimentally implemented 
as follows. We start with a spin-wave mode \a) with a 
real mode shape S(z), carry out backward retrieval, and 
measure the outgoing field. We then prepare the time 
reverse of the measured field shape and store it back into 
the ensemble using the time-reversed control pulse. The 
projections Pb and Pa happen automatically since we do 
not reverse the reservoir modes and the leak. The renor- 
malization can be achieved during the generation of the 
time-reversed field mode, while the time reversal for the 
spin wave will be unnecessary since a real spin wave will 
stay real under retrieval followed by time-reversed stor- 
age. The iteration suggested here is, thus, indeed equiv- 
alent to the iteration in Eq. (|2"T|) with 5(1 — z) replaced 
with S(z) (since Eq. (|2ip optimizes forward retrieval). 

For single photon states, the measurement of the out- 
going field involved in the procedure above will require 
many runs of the experiment. To circumvent this, one 
can use the fact that the equations of motion ©-(J?]) 
for the envelope of the quantum field mode are iden- 
tical to the equations of motion for the classical field 
propagating under the same conditions. One can, thus, 
use the optimization procedure with classical light pulses 
and find optimal pairs of input pulse shapes and control 
fields, which will give optimal storage into the spin wave 
Sd{l — z). However, since the equations of motion for 
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quantum light modes are identical to the classical prop- 
agation equations, this data can then be interpreted as 
optimal pairs of control fields and quantized input pho- 
ton modes for optimal storage of nonclassical light (such 
as single photons) into the optimal backward retrieval 
mode Sd(l — z). 

We will now briefly discuss the application of time re- 
versal ideas to the optimization of the combined process 
of storage followed by retrieval. For real spin waves, stor- 
age and backward retrieval are time reverses of each other 
since real spin waves are unaltered by complex conjuga- 
tion. Consequently, the time reversal iteration of storage 
and backward retrieval optimizes both of them, as well 
as the combined process of storage followed by backward 
retrieval. Therefore, for a given input, the storage con- 
trol field that optimizes storage alone will also be optimal 
for storage followed by backward retrieval. 

In contrast, (forward) storage and forward retrieval are 
not time reverses of each other, and the entire process of 
storage followed by forward retrieval has to be optimized 
as a whole. The general time reversal iteration procedure 
can still be used in this case with the understanding that 
the spaces A and B of initial and final states are the right 
propagating modes to the left of the ensemble (except 
for later empty input modes during retrieval) and right 
propagating modes to the right of the ensemble (except 
for earlier storage leak modes), respectively, while the 
remaining modes are reservoir modes, spin-wave modes, 
leak modes, and empty incoming photon modes from the 
left during retrieval. Since the time-reverse of storage 
followed by forward retrieval is itself storage followed by 
forward retrieval except in the opposite direction, the op- 
timization can be carried out physically by starting with 
a given input field mode, storing it and retrieving it for- 
ward with given control pulses, time reversing the output 
and the control pulses, and iterating the procedure. The 
optimal control-dependent input field, which the itera- 
tion will converge to, will then be stored into a particular 
optimal spin wave, which itself will be independent of the 
control used for the iteration. In Sees . I VI Cl and rVTTl we 
will look at storage followed by forward retrieval in more 
detail. 

It is important to note that the discussion in this sec- 
tion assumed that the two metastable states are degener- 
ate. If they are not degenerate, a momentum Afc = io sg jc 
will be written onto the spin wave during storage, so that 
its time reversal will no longer be trivial. In Sec. IVIII1 
we will discuss in detail how the optimization is modified 
when the metastable states are not degenerate. 

Procedures that are related to ours and that also use 
time-reversal iterations for optimization are a standard 
tool in applied optimal control [HI, [ill [2(| HH and 
have been used for a variety of applications in chem- 
istry [H, [H[, NMR [H, and atomic physics [H, HI]. 
In most of these works, time reversal iterations are used 
as a mathematical tool for computing, for a given ini- 
tial state | a), the optimal time-dependent control that 
would result in the final state with the largest projection 



on the desired subspace B of final states. In fact, this 
mathematical tool is directly applicable to our problem 
of shaping the control pulses, as we will discuss elsewhere 
[27l ]. However, our use of time reversal iterations in the 
present paper and in papers I and III differs in two ways 
from that of Refs. [IllialMllilSlMIMlHlIl- First, 
we use time reversal iterations to find the optimal \a) in 
the subspace A of initial states for a given propagator 
E/[fl(t)], rather than to shape the control Vl(t) itself for a 
given I a) (we shape the control by explicitly solving the 
equations, as explained in Sees . I VII and I VIip . Second, the 
time reversal iterations discussed in the references above 
are a purely mathematical technique, while our iterative 
algorithm can be implemented experimentally 

The main result of this section is an iterative procedure 
for solving or experimentally finding the optimal retrieval 
spin wave, while the main result of the previous section 
was that one can generate optimal pairs of inputs and 
control fields by time-reversing the output and the con- 
trol field of such optimal retrieval. In order to say, how- 
ever, for which input fields the optimal storage control 
O(t) can be found (or, equivalently, into which output 
fields can the optimal spin- wave excitation be retrieved), 
we need to consider the limits, in which Eq. (JT^J) can be 
fully solved analytically. These limits, adiabatic and fast, 
will be discussed in Sees. IVII and IVIll respectively. 

VI. ADIABATIC RETRIEVAL AND STORAGE 

A. Adiabatic retrieval 

1. Solution for the output field and for the retrieval 
efficiency 

Based on the branching ratio and the time reversal ar- 
guments, we have found the maximal storage efficiency 
at each d and have described the optimal storage sce- 
nario in the three preceding sections. Since a given input 
mode can be optimally stored if and only if optimal re- 
trieval can be directed into the time-reverse of this mode, 
in the following sections (Sees. I VII and IVII[) , we solve 
Eq. (|12|) analytically in two important limits to find out, 
which modes we can retrieve into and store optimally. 
The first of these two limits, which we will consider in 
the next five sections fSecs. lyTAl - IVIE)) . corresponds to 
smooth input and control fields, such that the term in the 
square brackets in Eq. (fT2| can be ignored. This "adi- 
abatic limit" corresponds to an adiabatic elimination of 
the optical polarization P in Eq. (fTTj) . The precise con- 
ditions for this adiabatic elimination will be discussed in 
Sec. IVIDI In this section (Sec. I VI A[) . we consider the 
retrieval process. 

Similar to the cavity discussion in Sec. V A of pa- 
per I, it is instructive to note that in the adiabatic ap- 
proximation (i.e., with dfP in Eq. ^ replaced with 0), 
rescaling variables £ and P by Q and changing variables 
t — » h{T r , t), where (as in Eq. (15) in paper I, except now 
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in dimensionless form) 



h(i, i') 



\n(t")\ 2 dt", 



(22) 



makes Eqs. (Jl])-© independent of f2. This allows one to 
solve these equations in an ^-independent form and then 
obtain the solution for any given f2 by simple rescaling. 
A special case of this observation has also been made 
in Ref. [28[, where the authors treat the Raman limit. 
However, since Eqs. ([5j)-([Zl) are relatively simple, we will 
avoid causing confusion by using new notation and will 
solve these equations without eliminating fl. 

To solve for the output field during adiabatic retrieval, 
we assume for simplicity that retrieval begins at time 
t = rather than at time t = T r and that the initial spin 
wave is S(z,t = 0) = S(z). In the adiabatic approxi- 
mation, Eqs. ([7]) and (fTl) reduce to a linear first order 
ordinary differential equation on S. Solving this equation 
for S(u,t) in terms of S(u'), expressing £(u,f) in terms 
of S(u, t) using Eqs. (fTT))) and (fTT|) . and taking the inverse 
Laplace transform u — > z = 1, we arrive at 



1 



£{l,t) = -VdQ(t) / dz 

Jo 1 + iA 



— e !+ iA 



.^aS si.-?. (23) 
V 1 + iA / 



The t-dependent and the z-dependent phases in the ex- 
ponent can be interpreted as the ac Stark shift and the 
change in the index of refraction, respectively. 
By using the identity [2!| 



drre-P r2 Io(Xr)Io^r) = I (^) (24) 



for appropriate [a, A, and p, we find that for a sufficiently 
large /i(0,oo) (dh(0,oo) 2> |rf + iA| 2 ), the retrieval effi- 
ciency (Eq. ^ with T r replaced with 0) is 

Vt= f dz f dz'k T (z,z')S(l -z)S*(l-z'), (25) 
Jo Jo 

in agreement with Eq. (fbl)) . So r) r is independent of de- 
tuning and control pulse shape but depends on the spin 
wave and the optical depth. Thus, the adiabatic approx- 
imation does not change the exact value of efficiency and 
keeps it independent of detuning and classical control 
field by preserving the branching ratio between desired 
and undesired state transfers. It is also worth noting 
that, in contrast to Eq. (|14|) . Eq. (j23|) allows us to treat 
and optimize retrieval even when the energy in the con- 
trol pulse is limited [1| (i.e., dh(0,oo) < \d + iA\ 2 ). 
However, in the present paper, the treatment of adia- 
batic retrieval is focused on the case when the control 
pulse energy is sufficiently large (dh(0, oo) 3> \d + iA\ 2 ) 
to leave no excitations in the atoms and to ensure the 
validity of Eq. (or, equivalently, Eq. (fTI)) ). 



As pointed out in the Introduction to paper I, two im- 
portant subsets of the adiabatic limit, the resonant limit 
and the Raman limit, are often considered in the litera- 
ture because the equations can be simplified in these lim- 
its. Although we demonstrate in this work that the basic 
underlying physics and hence the optimal performance 
are the same for these two photon storage techniques, a 
more detailed analysis reveals significant differences. It 
is precisely these differences that obstruct the underlying 
equivalence between the two protocols. And it is these 
differences that make this equivalence remarkable. As an 
example of such a difference, resonant and Raman limits 
give different dependence on d of the duration T out of the 
output pulse. To see this, it is convenient to ignore the 
decay in Eq. (|23|) (due to the rescaling, this means we 
ignore 1 in 1 + iA). If we do this, we obtain 



£(l,t) = iVdhji) J o dzie iM2 ^ Jq U ^ )d ~ Z - 



xS(l-z), 



(26) 



where Jq(x) = Io(ix) is the zeroth order Bessel function 
of the first kind. In the resonant limit (c^y 3> |A|), we 
find the duration of the output pulse by observing that 
the A -> limit of Eq. ([26]) is 



~V1 



s 1 



h(0, t) 



(27) 



with the understanding that S(z) vanishes outside of 
[0,1]. This is just the ideal lossless and dispersionless 
group velocity propagation of Refs. [H, H3|, also know 
as EIT. This implies a duration T out ~ rf'y / 1 1 2 for the 
output pulse in the resonant limit, which is consistent 
with the cavity case analyzed in paper I if one identifies 
C and d. In the Raman limit (d*y <C |A|), the length of 
the output pulse is simply given by the fall-off of Jo and 
is found from h(Q, i)d/A 2 - 1 to be T out ~ A 2 /( 7 d|f7| 2 ). 
It is worth noting that the appropriate Raman limit con- 
dition is "fd <C |A| and not 7 <C |A| as one might naively 
assume by analogy with the single-atom case. It is also 
important to note that if one is limited by laser power 
(as in Ref. [28|) and desires to achieve the smallest pos- 
sible T out , the above formulas for T out imply that EIT 
retrieval is preferable over Raman. 



2. Dependence of retrieval error on optical depth d 

In the cavity case analyzed in paper I, only one spin- 
wave mode is available and the retrieval error is always 
1/(1 + C) (ssl/C for C ^> I). In free space, in contrast, 
the retrieval error depends on the spin-wave mode S(z) 
and as we will explain in this section, scales differently 
with d depending on the spin wave. Since the retrieval 
efficiency is independent of A, to gain some physical in- 
tuition for the error dependence on the spin wave and on 
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d, we will focus on the A = case, for which the for- 
malism of EIT transparency window has been developed 
[3lj ] . For A = and large d, the integrand in Eq. |23|) 
can be approximated with a Gaussian. Then using di- 
mensionless momentum k — kL and defining the Fourier 
transform of S(z) as S(k) = (2n)^ 1 dzS(z) exp(— ifcz), 
we can write Eq. (|23|) as 



£{1,1) 



h(t) 

~7d J- 



dke Ml 



(28) 

In the limit d — * oo, the Gaussian term can be replaced 
with 1 to yield back the group velocity propagation in 
Eq. (j27|) . Computing the efficiency using Eq. (|28]) , we 
find, after a change of variables t — > r = h(0, t)/d, 



Vr 



dr 



dk/ k(1 - T) e'^S{k) 



(29) 



We will now show that the Gaussian term of width 
A/ceit = \J~dpr in the integrand in Eq. (|29| can be in- 
terpreted as the effective momentum-space EIT trans- 
parency window for the spin wave. We start by noting 
that the equivalent of t in the original units (call it z prop ) 
is equal to z pr0 p = Lr{t) = J Q v g {t')dt' and, thus, repre- 
sents the propagation distance (v g — cQ 2 /(g 2 N) is the 
EIT group velocity (3lj). This interpretation of z prop 
also follows from the fact that, if one ignores the Gaus- 
sian in Eq. (|2"9"j). the spin wave would be evaluated at 
z = 1 — t. Thus, in terms of the propagation distance 
•Zprop, the width of the momentum-space transparency 
window in Eq. (|29p can be written, in original units, 
as Afc E iT = Ak ElT /L = yfdfr/L = ^/ g 2 N/ (~/cz pvop ). 
Thus, as the propagation distance z plop decreases, the 
width A/ceit of the transparency window gets wider and 
eventually becomes infinite at the z = 1 end of the ensem- 
ble, where z prop = 0. The consistency of our expression 
AfcEiT for the effective momentum-space EIT window 
with the expression for the frequency-space EIT trans- 
parency window Aojeit = v g \J g 2 N / (7Cz prop ) [3l| imme- 
diately follows from rescaling by v g both Aweit and the 
dark state polariton [12t |30j bandwidth Aw p = v g Afc sp ; n 
(where Afc sp ; n is the width of S(k) in the original units). 
In fact, the change of variables t — ► r that led to Eq. (j2"9"|) 
precisely accomplished this rescaling of the polariton and 
the EIT window by the group velocity. It is worth not- 
ing that this proportionality of both the polariton band- 
width and the frequency-space EIT window width to the 
group velocity (and, hence, the existence of the control- 
independent effective momentum-space EIT window) is 
another physical argument for the independence of re- 
trieval efficiency from the control power. 

An important characterization of the performance of 
an ensemble-based memory is the scaling of error with op- 
tical depth at large optical depth. In the cavity case ana- 
lyzed in paper I, there was only one spin- wave mode avail- 
able and the retrieval error for it was 1/(1 + C) (« 1/C 
for C> 1), where the cooperativity parameter C can be 



thought of as the effective optical depth enhanced by the 
cavity. By qualitative arguments, we showed in Sec. IIIII 
that the retrieval efficiency in free space is 1 — rj ~ 1 jd. A 
more precise value can be found numerically from the op- 
timal spin wave fSec. lIII[ ). which gives a maximal retrieval 
efficiency that scales approximately as ~ 2.9/d, i.e., one 
over the first power of density, precisely as in the cavity 
case. However, this 1/d scaling turns out to be not the 
only possibility in free space. The scaling of the retrieval 
error with d can be either 1/d or 1/yd depending on the 
presence of steps (i.e., discontinuities in the amplitude or 
phase of the spin wave). Specifically, numerics show that 
for a spin wave that does not have steps at any z < 1, 
the retrieval error scales as 1/d, while steps in the phase 
or amplitude of the spin wave result in a 1/ \fd error. In 
particular, a step in the amplitude of S(z) at position z of 
height I can be found numerically to contribute an error 
of I 2 \j2/-K\/1 — z/Vd at large d. The reason for the im- 
portance of steps is that at high d, the effective EIT win- 
dow is very wide and only the tails of the Fourier trans- 
form S(k) of the spin wave S(z) matter. A step, i.e., a 
discontinuity, in the function S(z) means that its Fourier 
transform falls off as S(k) ~ 1/k. Thus, if we assume all 
frequencies outside of an effective EIT window of width 
A^eit = \fdpr ~ y/d/(l — z) get absorbed, the error 
will be proportional to dk\S(k)\ 2 ~ \/ (1 — z)/d, 

precisely as found with numerics. Numerics also show 
that if a step in |<S(5)| is not infinitely sharp, at a given 
d, a feature should be regarded as a step if the slope of 
^(z)! 2 is bigger than ~ \/d (and will then contribute 
a 1/ \fd error) . The simple physical reason for this is 
that only if a feature in S(z) is narrower than 1/yd will 
it extend in fc space outside the effective EIT window of 
width A^eit = ydjr. While we only performed detailed 
analysis of steps in the amplitude of S(z), steps in the 
phase of S(z), as we have already noted, also contribute 
a l/Va error, and we expect that similar dependence on 
the position and sharpness of such phase steps holds. 

A useful analytical result on scaling that supports these 
numerical calculations is the error on retrieval from a flat 
spin wave S(z) = 1, which can be calculated exactly from 
Eq. (25JI to be 



l-Vr = e- d {I (d) + h(d)). 



(30) 



Using the properties of modified Bessel functions of the 
fir st ki nd, one finds that as d — > oo, the error approaches 
y/2/w/y/d, which is consistent with the general formula 
since a fiat spin wave has one step at z < 1, i.e., a step 
of height 1 at z = 0. In fact, it is this analytical result 
that allows us to exactly identify the \Z2/tc prefactor in 
the error due to amplitude steps. 

Based on the results of this section and the results of 
paper I, we can identify several advantages of using a 
cavity setup. First, in the cavity, the optical depth is 
enhanced by the value of the cavity finesse from the free- 
space value of d to form the cooperativity parameter C . 
Moreover, in terms of d and C the errors during opti- 
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mal storage in free space and in the cavity scale as 2.9/d 
and 1/C, respectively. That is, even if one ignores the 
enhancement due to cavity finesse, the cavity offers a fac- 
tor of 3 improvement. In addition to that, if one is forced 
to retrieve from the flat spin-wave mode S(z) = 1 (which 
is the case, for example, if the spin wave is generated via 
spontaneous Raman scattering as in Ref. [5J), the free- 
space error is increased from the optimal and scales as 
y/2/ir/ \fd 1 while in the cavity case the mode S(z) = 1 
is, in fact, the only mode coupled to the cavity mode and 
is, therefore, precisely the one that gives 1/C scaling. 
Finally, because there is only one spin wave mode acces- 
sible in the cavity setup, the time reversal based iterative 
optimization procedure (Sec. EJ requires only one itera- 
tion in the cavity case. On the other hand, the free-space 
setup described in this paper is much simpler to realize 
in practice and allows for the storage of multiple pulses 
in the same ensemble, e.g., time- bin encoded qubits [32| . 

Using the effective EIT window concept developed in 
this section, we can now interpret from a different per- 
spective the optimal retrieval spin waves computed in 
Sec. IIIII These spin waves represent at each d the opti- 
mal balance between maximal smoothness (to minimize 
the momentum space width Afc sp i n of the spin wave so 
that it better fits inside the effective EIT window) and 
least amount of propagation (to minimize r and, thus, 
maximize the width A/ceit = ydfr of the effective EIT 
window itself). 



3. Shaping retrieval into an arbitrary mode 

We have shown that optimal storage of a given input 
mode requires the ability to retrieve optimally into the 
time reverse of this input mode. Thus, by finding the 
modes we can retrieve into, we will also find the modes 
that can be optimally stored. In this section we prove 
that by adjusting the control during retrieval we can re- 
trieve from any mode S[z) into any given normalized 
mode £2(^)1 provided the mode is sufficiently smooth to 
satisfy the adiabaticity condition (which in the original 
units means T out <i7 3> 1, where T QUt is the duration of 
£2^), as we discuss in Sec. IVID|) [33j j . 

We know from Sec. IIIII that the retrieval efficiency rj T 
is independent of the detuning A and the control CI, pro- 
vided the retrieval is complete (for adiabatic retrieval, 
the condition on the control pulse energy for complete re- 
trieval is dh(0,oo) > \d + iA| 2 , as found in Sec. lVIAl) . 
Thus, to find the control that retrieves S(z) into any 
given normalized mode E-zyt) with detuning A, we need 
to solve for (lit) in equation (j23[) with £(1, t) — y/fji£-2.(t)'. 



ifr£a(*) = -Vdn(t) / dz 

1 -MA 



I _ h(0,t) + dz 



x/o 2«S S(l-i). (31) 
V 1 + iA / 



To solve for Cl(t), we integrate the norm squared of both 
sides from to t and change the integration variable from 
r to h! = h(0,t') on the right hand side to obtain 



Vl 1 dt'\£ 2 (t')\ 2 = 
lo 



h(0,t) 



dti 



1 , Vd h '+ d * 

dz -e !+ iA 

1 + iA 

2 



x/ 



. 2 rS *'-*> 



(32) 



In the cavity case, the corresponding equation (20) of 
paper I was solvable analytically. To solve Eq. ([32)) for 
h(0, t) numerically, we note that both sides of Eq. (f3"2"|) are 
monotonically increasing functions of t and are equal at 
t = and t = 00 (provided /i(0, 00) can be replaced with 
00, which is the case if dh(0, 00) 3> |<i + iA| 2 ). Therefore, 
Eq. (|32p can always be solved for h(0,i). \Cl(t)\ is then 
deduced by taking the square root of the derivative of 
h(0,t). The phase of fl is found by inserting |f2| into 
Eq. (|3"Tj) and is given by 



Arg 



7r + Arg \£ 2 ( 



1 + A 2 



-Arg 



dz— 
1 + iA 



— e !+ iA 



The second and third terms are the phase of the desired 
output and the compensation for the Stark shift, respec- 
tively. In the resonant limit (d*y ^S> |A|), we can set 
A = 0. Then assuming the phase of S(z) is indepen- 
dent of z, the phase of the optimal control is given (up 
to a constant) solely by the phase of the desired output. 
In the Raman limit [dry -C |A|), the integral in the last 
term is approximately real but at times near the end of 
the output pulse (t « T out ) it can change sign and go 
through zero. At these times, the optimal \Cl(t)\ diverges 
and the phase of O(t) changes by tt. Numerical simu- 
lations show, however, that can be truncated at 
those points without significant loss in the retrieval effi- 
ciency and in the precision of £-i(t) generation. Moreover, 
these points happen only near the back end of the de- 
sired output pulse over a rather short time interval com- 
pared to the duration of the desired output pulse. We 
can therefore often even completely turn off Cl(t) during 
this short interval without significantly affecting the re- 
sult, so that the problem of generating large power and tt 
phase shifts can be avoided altogether (see Sec. IVIDI for 
another brief discussion of this issue). However, these po- 
tential difficulties in the Raman limit for generating the 
optimal control (which also has to be chirped according 
to Eq. (|3"3"|) in order to compensate for the Stark shift) 
make the resonant (EIT) limit possibly more appealing 
than the far-off-resonant (Raman) limit. 

Finally, we note that a divergence in |0(t)| can occur 
at any detuning A even when the above Raman-limit 
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divergences are not present. Specifically, similarly to the 
cavity discussion in paper I, if one wants to shape the 
retrieval into a mode £2\t) that drops to zero at some 
time T out sufficiently rapidly, |f2(t)| will go to oo at t = 
T ou t. However, as in the above case of the Raman- limit 
divergences, the infinite part can be truncated without 
significantly affecting the efficiency or the precision of 
82(b) generation. One can confirm this by inserting into 
the adiabatic solution in Eq. ([23]) a control pulse that is 
truncated to have a value of /i(0, 00) that is finite but 
large enough to satisfy dh(0, 00) ^> \d + iA\ 2 . However, 
to be completely certain that the truncation is harmless, 
one has to solve Eqs. ©-([7]) numerically without making 
the adiabatic approximation. We will do this in Sec. lVI Dl 
for the case of storage, where the same truncation issue 
is present. 



B. Adiabatic storage 

In principle, the retrieval results of the previous sec- 
tion and the time reversal argument immediately imply 
that in the adiabatic limit (see Sec. lVIDl for precise con- 
ditions), any input mode £i n (t) at any detuning A can be 
stored with the same ci-dependent maximum efficiency if 
one appropriately shapes the control field. However, for 
completeness and to gain extra physical insight, in this 
section, we present an independent solution to adiabatic 
storage. 

Using the Laplace transform in space and a procedure 
similar to the one used in Sec. lVI Al to solve retrieval, we 
find that the adiabatic solution of storage is 

1— f T ~~ 1 h(i,f) + dz 

S(z,T) = -Vd / dtSl*{t) i+*a 

Jo 1+iA 



Jh(t,T)dz 

x/o I 2 V =— I & m (t). (34) 

1 + iA 



It is important to note that the retrieval equation (|23|) 
and the storage equation (f34|) can be cast in terms of the 
same ^-dependent function m as 



£ out (i) 



dzm[n(t'),t,z] 5(1-5), (35) 



S(z,f) = / dtmQ*(T -t'),T -t,z £„(£). (36) 

This is precisely the general time reversal property of our 
equations that we discussed abstractly in Sees. HVI andfVl 
and used to find the optimal storage strategy from op- 
timal retrieval. However, as we said in the beginning of 
this section, we will now, for completeness, optimize stor- 
age directly without using time reversal and our solution 
for optimal retrieval. 

We would like to solve the following problem: given 
£in{t), A, and d, we are interested in finding Cl(t) that 



will give the maximum storage efficiency. To proceed 
towards this goal, we note that if we ignore decay 7 and 
allow the spin wave to extend beyond z = 1, we get the 
"decayless" storage equation 



s(z) = / dtq(z,i)S- m (T), 



(37) 



where the "decayless" mode s(z) is defined for z from 
to 00 instead of from to 1 and where 



q{z,i) = i%/dn*(i)— e 4 



Jo 2 



h(t, T)dz 

"a 



(38) 

Since in Eq. l[3"7|). both sources of storage loss (the de- 
cay 7 and the leakage past z = 1) are eliminated, the 
transformation between input modes £i n (t) and decay- 
less modes s(z) becomes unitary. Indeed, we show in 
Appendix [E] that Eq. (I37|) establishes, for a given f2(t), 
a 1-to-l correspondence between input modes £i n (t) and 
decayless modes s(z). Moreover, we show in Appendix 
[E]that Eq. ([37]) also establishes for a given s(z) a 1-to-l 
correspondence between input modes £ m (t) and control 
fields Cl(t). In particular, this means that we can com- 
pute the control field that realizes decayless storage (via 
Eq. ([3"T]) ) of any given input mode £ ln (t) into any given 
decayless spin- wave mode s(z). 

A key element in the control shaping procedure just 
described is the ability to reduce the problem to the uni- 
tary mapping by considering the decayless (and leakless) 
solution. The reason why this shaping is useful and why 
it, in fact, allows us to solve the actual shaping prob- 
lem in the presence of decay and leakage is that the spin 
wave, into which we store in the presence of decay, can 
be directly determined from the decayless solution: us- 
ing Eqs. ([Ml) and ([27]) and Eq. (J21J (with appropriate 
complex values of /i, A, and p), we find that 



S(z,f) 



dz'de- d{i+i ' ) I (2dVl~z 1 )s{~z'). (39) 



This means that, remarkably, S(z, T) (and hence the 
storage efficiency) depends on £- m (t) and Q(t) only 
through the decayless mode s(z), which itself can be com- 
puted via unitary evolution in Eq. (j37[) . 

Computing storage efficiency from Eq. ([39]) as a func- 
tional of s(z) and maximizing it under the constraint 
that s(z) is normalized gives an eigenvalue problem sim- 
ilar to Eq. p6p except the upper limit of integration is 
00 and the kernel is different. After finding the optimal 
s(z) via the iteration scheme similar to the one used to 
solve Eq. I|16p . we conclude the procedure for optimal 
storage control shaping by using the unitary transforma- 
tion in Eq. ([37| to solve for the control in terms of s(z) 
and £- m (i), as shown in Appendix [El Since efficiency is 
determined by s(z) alone, this gives the optimal stor- 
age with the same maximal efficiency for any input pulse 



shape (provided it is sufficiently smooth, as discussed in 
Scc. rvTDl) . 

Having derived the optimal storage control, we can now 
explicitly verify the results obtained from the time rever- 
sal reasoning. In Appendix [El we show that the mode 
S(z, T) used in optimal storage is just the optimal mode 
for backward retrieval; that the optimal storage efficiency 
and optimal retrieval efficiency are equal; and that the 
optimal storage control for a given input mode is the 
time-reverse of the control that gives optimal backward 
retrieval into the time-reverse of that input mode. 

To give an example of optimal controls, we consider a 
Gaussian-like input mode (shown in Fig. [3]) 

£ in (t) = A{e^/ f -^ 2 - e- J - 5 )/Vf, (40) 

where for computational convenience we set £i n (0) = 
£in{T) = and where A w 2.09 is a normalization con- 
stant. Figure [3] shows the corresponding optimal storage 
control shapes for the case A = and d = 1, 10, 100, 
as well as the limiting shape of the optimal Q as d — ► oo. 
The controls are plotted in rescaled units so that the 
area under the square of the curves shown is equal to 
J T dt'v g (t') (where v g is the EIT group velocity), 
which is the position (in units of L), at which the front 
end of the pulse would get stored under ideal decayless 
propagation. From time reversal and the condition for 
complete retrieval, it follows that the control pulse en- 
ergy (oc h(0,T)) and hence L^ 1 dt'v g (t') should di- 
verge. Thus, at any finite d, the optimal f2 plotted in 
Fig. [3] should actually diverge at t = 0. However, the 
front part of the control pulse affects a negligible piece 
of £in(i), so truncating this part (by truncating s(z) at 
some z, for example) does not affect the efficiency (we 
will confirm this again in Sec. |yi"D|) , Naively, the opti- 
mal control should roughly satisfy L _1 J Q dt'v s (t') = 1: 
the control (and, thus, the dark state polariton group 
velocity) should be small enough to avoid excessive leak- 
age; and at the same time it should be as large as possible 
to have the widest possible EIT transparency window to 
minimize spontaneous emission losses. For the truncated 
optimal controls, we see that L~ x J Q T dt'v g (t') (i.e., the 
area under the square of the curves in Fig. [3]) is, in fact, 
greater than 1. As d decreases, L~ x J Q dt'v s (t') decreases 
as well, and allows for less and less leakage so that only 
as d — > oo, s(z) — > \/3(l — z), L^ 1 J Q T dt'v g (t') — > 1, and 
no leakage is allowed for. 

Similarly to the cavity storage discussed in paper I, 
although optimal storage efficiencies are the same in the 
Raman and adiabatic limits, the two limits exhibit rather 
different physical behavior. It is now the dependence on 
d of the intensity of the optimal control field that can be 
used to distinguish between the resonant and the Raman 
regimes. Using an analysis very similar to pulse duration 
analysis of Sec. IVI Al or alternatively relying on the fact 
that optimal storage is the time-reverse of optimal re- 
trieval, we find that in the resonant limit (dj S> |A|), 
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FIG. 3: Input mode £in(t) (dashed) defined in Eq. (|40|l and 
control fields Q(t) (in units of \J d-y/T) that maximize for this 
£in(t) the efficiency of resonant adiabatic storage (alone or 
followed by backward retrieval) at d — 1, 10, 100, and d — > oo. 



|£!| 2 ~ drf/T, while in the Raman limit (dry < |A|), 
|£!| 2 ~ A 2 /(jTd). Both of these agree with the cavity 
model of paper I (identifying C and d) while the resonant 
control agrees with v g T ~ L. As explained in the cavity 
case in Sec. V B of paper I, this opposite dependence 
of |f2| on d in the Raman and EIT limits is the conse- 
quence of the fact that the coupling of the input photon 
to the spin wave has the opposite dependence on £1 in the 
two regimes. Finally, we note that in the Raman limit 
(|A| > dry), \n\ 2 ~ A 2 /(yTd) > dy/T, which means 
that if one is limited by control power (as in Ref. [28|), 
the EIT regime is preferable to the Raman regime. 



C. Storage followed by retrieval 

In the cavity case discussed in paper I, there was only 
one spin-wave mode we could write on. Moreover, this 
spatially uniform mode looked the same in the forward 
and backward directions (assuming negligible metastable 
state splitting u) sg ). Therefore, optimal storage into that 
mode guaranteed that the combined procedure of storage 
followed by retrieval was optimal as well and the total ef- 
ficiency did not depend on the retrieval direction. In 
contrast, the free-space model allows for storage into a 
variety of modes, each of which has a different retrieval 
efficiency that is also dependent on retrieval direction. 
Therefore, in free space, we will first discuss the opti- 
mization of storage followed by backward retrieval and 
then the optimization of storage followed by forward re- 
trieval. 

Since we have shown that the optimal spin- wave mode 
for backward retrieval is also the optimal mode for stor- 
age, the controls found in Sec. IVI Bl which optimize stor- 
age, are also optimal for storage followed by backward 
retrieval. Figure 0] shows the maximum total efficiency 
for storage followed by backward retrieval {rj^^ - solid 
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FIG. 4: 7?bSS (solid) and r/J-^ (dotted) are maximum total 
efficiency for adiabatic (or fast) storage followed by backward 
or forward retrieval, respectively. Square (dashed) is the total 
efficiency for resonant storage of £i u (t) from Eq. ((40} followed 
by backward retrieval, where the storage control field is a 
naive square pulse with the power set by v s T — L. 



line), which in the adiabatic limit can be achieved for 
any input pulse. For comparison, we also show the total 
efficiency for storage followed by backward retrieval for a 
Gaussian-like input mode defined in Eq. (|40p (assuming 
the adiabatic limit Td"f 3> 1) with naive square stor- 
age control pulses on [0, T] with power set by v g T = L, 
where u g is the EIT group velocity (?7squarc - dashed line) . 
The significant increase in the efficiency up to the input- 
independent optimal efficiency r^ck due to the use of 
optimal storage control pulses instead of naive ones is, 
of course, not unique to the input pulse of Eq. (|40[) and 
holds for any input pulse. In fact, since at moderate 
values of d the naive control pulse is far from satisfying 
the complete retrieval condition, it is not optimal for any 
input. 

Since the optimal mode for storage or retrieval alone 
is not symmetric, a separate optimization problem has 
to be solved for the case of storage followed by forward 
retrieval. We show in Appendix [E] that Eq. (|57|) sets 
up, for any sufficiently smooth £- m (t), a 1-to-l correspon- 
dence between decayless modes s(z) and control fields 
n(t). Moreover, the decayless mode alone determines the 
total efficiency of storage followed by forward retrieval, 
which can be found by inserting Eq. ([39]) into Eq. (|25)l . 
Thus, the optimization problem reduces to finding the 
optimal s(z) by the iterative optimization procedure ex- 
cept with a new kernel. 

From a different perspective, since the process of stor- 
age followed by forward retrieval as a whole fits into the 
general setup of Sec. |V] we can use time reversal to op- 
timize it. In particular, we showed in Sec. fVlhow to find 
the optimal input mode to be used with a given pair of 
control pulses: we should take a trial input mode, carry 
out storage followed by forward retrieval, time reverse 



the output and the controls, and iterate until conver- 
gence. Since the reverse process is itself storage followed 
by forward retrieval (except in the opposite direction) 
and since the iterations will optimize it as well, the spin- 
wave mode used in optimal storage followed by forward 
retrieval must be the one that flips under forward re- 
trieval, followed by time reversal and (backward) stor- 
age. Moreover, it follows that the control pulse that we 
should use for a given £{ n (t) is the time- reverse of the 
control that retrieves the flipped version of the optimal 
spin- wave mode backwards into £* n {T — t). 

Thus, instead of computing the optimal s(z), we can 
solve the following eigenvalue problem that finds the op- 
timal mode to store into 



XS{z) 



dz'k T (z, l-z')S(z'). 



(41) 



This eigenvalue equation is just a simple modification of 
the retrieval eigenvalue equation (|16p : we are now com- 
puting the mode that flips under forward retrieval fol- 
lowed by time reversal and storage, while Eq. (fTS)) finds 
the mode that stays the same under backward retrieval 
followed by time reversal and storage. The total effi- 
ciency of storage followed by forward retrieval is then A 2 . 
However, in contrast to storage followed by backward re- 
trieval, the storage efficiency and the forward retrieval ef- 
ficiency during the optimal procedure are not each equal 
to A; the storage efficiency is greater. 

The optimal spin-wave modes that result from Eq. (|4"Tj) 
are shown in Fig. [5] for the indicated values of d. At 
small d, the optimal mode is almost flat since at small d 
the optimal retrieval and storage modes are almost flat 
and, thus, almost symmetric. As d increases, the opti- 
mal mode first bends towards the wedge -\/3(l — z) simi- 
larly to the optimal backward retrieval mode. But then 
above d ft 3, the optimal mode starts shaping towards 
the parabola S{z) = ^15/8(1 - 4(z - 1/2) 2 ), which, as 
expected, avoids the 1/Vd error from discontinuities by 
vanishing at the edges and, simultaneously, maximizes 
smoothness. 

The maximal total efficiency for storage followed by 
forward retrieval fy^w ' s plotted as a dotted curve in 
Fig- SI C (~ 1 - 19/d as d — > oo) is less than ry™^ 
1 — 5.8/d as d — * oo) since for optimal backward re- 
trieval, storage and retrieval are each separately optimal, 
while for forward retrieval a compromise has to be made. 
From a different perspective, forward retrieval makes it 
more difficult to minimize propagation since the whole 
excitation has to propagate through the entire medium. 



D. Adiabaticity conditions 

We have shown that in the adiabatic limit, any input 
mode can be stored with the same maximum efficiency. 
In this section, we show that independent of A, the neces- 
sary and sufficient condition for optimal adiabatic stor- 
age of a pulse of duration T to be consistent with the 
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FIG. 5: For different values of d, the optimal spin-wave mode 
to be used during storage followed by forward retrieval. 



adiabatic approximation is 

Td 1 > 1, 



(42) 



which is identical to the corresponding condition in the 
cavity model in paper I, except with C replaced with d. 
In fact, we omit the derivation of this condition since the 
argument goes along the same lines as the correspond- 
ing argument in the cavity case in Sec. V C of paper I, 
provided one uses the fact that in Eq. (fTTjl of the present 
paper |u| ~ 1 (since u is the Laplace variable correspond- 
ing to z, which, in turn, runs from to 1). 

Therefore, we immediately turn to the tasks of verify- 
ing numerically that Eq. (14"2")) is indeed the correct adia- 
baticity condition and of investigating the breakdown of 
adiabaticity for short input pulses. We consider adiabatic 
storage of a Gaussian-like input mode defined in Eq. (140]) 
and shown in Fig. [3l We use our adiabatic equations 
to shape the control pulse but then compute the total 
efficiency of storage followed by backward retrieval nu- 
merically from Eqs. (J5|)-([7|) without making the adiabatic 
approximation. As Tdj decreases to 1 and below, we ex- 
pect the efficiency to be reduced from its optimal value. 
In Fig. [6](a), the total efficiency of this procedure is plot- 
ted as a function of Td^ for A = and d = 1, 10, 100, 
and 1000. The dashed lines are the true optimal effi- 
ciencies. As expected, when Td^j < 10, the efficiency 
drops. In Fig.[6](b), we fix d = 10 and show how optimal 
adiabatic storage breaks down at different detunings A 
from to 2OO7. As in the cavity case of paper I, we see 
from Fig. &[b) that as we move from the resonant limit 
(dy 3> |A|) to the Raman limit (dj <C |A|), we can go to 
slightly smaller values ofTdj before storage breaks down. 
However, since the curves for A = IOO7 and A = 2OO7 
almost coincide, Td^ ^ 1 is still the relevant condition 
no matter how high A is. 

Before concluding the discussion of adiabaticity con- 
ditions, we note that, exactly as in the cavity case dis- 
cussion in Sec. V C of paper I, the magnitudes |f2(f)| of 
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FIG. 6: Breakdown of optimal adiabatic storage in free space 
at Td'y < 10. In (a), the total efficiency of storage followed 
by backward retrieval is plotted for A = and d = 1, 10, 100, 
and 1000. The horizontal dashed lines are the maximal values. 
Dotted lines are obtained for the input from Fig. [3] using adi- 
abatic equations to shape the storage control but then using 
exact equations to numerically compute the efficiency. In (b) , 
the same plot is made for d = 10 and A/7 = 0, 1, 10, 100, 200. 



the optimal control pulses, which go to 00 at t = 0, had 
to be trucated to generate Fig. [6] Moreover, Raman- 
limit divergences discussed (for the case of retrieval) 
in Sec. IVI A 31 had to be truncated at two places near 
t = 0.2 T for A/7 = 100 and 200 in Fig. |6£b). As in Sec. 
V C of paper I, the fact that the optimal efficiency given 
by the dashed lines in Fig. [5] is achieved by the dotted 
curves (obtained with truncated controls) is the proof 
that truncation of the storage controls does not signif- 
icantly affect the storage efficiency. Since time reversal 
discussed Sec. II VI relates retrieval to optimal storage, this 
also means that truncation of the retrieval control fields 
does not significantly affect the precision, with which a 
given output mode £2$) can be generated. The losses 
associated with truncation are insignificant only if condi- 
tion dh(0, 00) > \d + iA\ 2 (see Sec. IVI A lj) is satisfied for 
the truncated retrieval control pulse (the same condition 
with 00 replaced with T applies to storage) . If the energy 
in the control pulse is so tightly limited that this condi- 
tion is violated, a separate optimization problem has to 
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be solved. This problem has been considered in Ref. [28| 
for the special case of Raman storage in free space in the 
limit of negligible spontaneous emission loss. Although 
Eqs. (f25| and (j3~4"|) allow one to treat and optimize the 
case of limited control pulse energy, this problem is be- 
yond the scope of the present paper. 



E. Effects of nonzero spin-wave decay 

We have so far assumed that the spin- wave decay rate 
in Eq. © is negligible (i.e., 7 S = 0). If this decay rate 
is not negligible but the corresponding incoming noise 
is vacuum (which is often true experimentally, as noted 
in Appendix [A] of the present paper and as explained in 
detail in Appendix A of paper I), the effect of nonzero 
spin- wave decay can be included simply by adding a term 
— %S to the right-hand side of Eq. (|7|) (where 7 S = 73/7). 
In this section, we discuss the effects such a term has on 
the adiabatic solution and optimal control field shaping 
discussed in Secs. lVI AllVIBI and I VI CI We will show, in 
particular, that nonzero 7 S simply introduces exponen- 
tial decay into the retrieval and storage solutions with- 
out making the optimization harder. We will also show 
that with nonzero spin-wave decay, optimal efficiencies 
become dependent on the input mode (or, equivalently, 
on the control fields). 

We first consider adiabatic retrieval discussed in 
Sec. IVI Al One can easily check that nonzero spin- 
wave decay rate 7 S simply introduces decay described by 
cxp(— %t) into Eq. (j!?3"|) , and, unless we retrieve much 
faster than 1/7 S , this makes the retrieval efficiency con- 
trol dependent. Moreover, if a given fixed control field is 
not strong enough to accomplish retrieval in a time much 
faster than 1/7 S , the problem of finding the optimal re- 
trieval mode for this particular retrieval control will give 
a different answer from the 7 S = case. In particular, 
for forward retrieval, as we increase 7 S (or, alternatively, 
decrease the power of the retrieval control), to minimize 
propagation time at the cost of sacrificing smoothness, 
the optimal retrieval mode will be more and more con- 
centrated towards the z = L end of the ensemble. As in 
the 7 S = case, we can find these optimal modes either 
by computing the (now ^-dependent) kernel to replace 
k T in Eq. (f!?5")) and its eigenvector with the largest eigen- 
value or, equivalently, by doing the iteration of retrieval, 
time reversal, and storage. 

The inclusion of nonzero 7 S also does not prevent us 
from being able to shape retrieval to go into any mode, as 
described for 7 S = in Sec. IVI A 31 We should just shape 
the control according to Eq. (|31|) as if there were no spin 
wave decay except the desired output mode on the 
left hand side should be replaced with the normalized 
version of £2^) exp(7 s i), i.e., 



The retrieval efficiency will, however, be output-mode- 
dependent in this case: it will be multiplied (and hence 

reduced) by a factor of [/ °° di'\£ 2 {t')\ 2 exp^i')] ~\ 
Since this factor is independent of the spin wave, even 
with nonzero 7 S , the optimal retrieval into £2{t) is 
achieved by retrieving from Sd(z) . 

We now turn to adiabatic storage discussed in 
Sec. IVI Bl One can easily check that nonzero spin-wave 
decay 7 S simply introduces exp(— %(T — t)) decay into 
Eq. (|38p (or, equivalently, into the integrand on the right 
hand side of Eq. J34])). Eq. (|39|) holds even with nonzero 
7 S . The optimal storage control can then be found from 
Eq. (|37[) as if there were no decay but the input mode 
were replaced according to 



£in(t) -> £in{t)e 



-%(T-t) 



dt'\£ m {t')\ 2 e 



2 -2%(f-i') 



(44) 

However, the optimal storage efficiency will now depend 
on input pulse duration and shape: it will be multiplied 

(and hence reduced) by dt'\£ la (t')\ 2 exp(-2%(f-t')). 
It is also important to note that nonzero 7 S still keeps the 
general time reversal relationship between storage and 
retrieval exhibited in Eqs. (|35|) and (|36jl . However, for 
£ 2 (t)=£* n (f-t), 



dt\£ in (t)\ 2 e- 2 ^ f -^ 



> 



dt\£ 2 (t)\ z e 



2„27 s t 



£ 2 (t)^£ 2 (t)e 



dt'\£ 2 (t' 



(43) 



(45) 

which means that with nonzero 7 S , the optimal storage ef- 
ficiency of a given input mode is greater than the optimal 
retrieval efficiency into the time-reverse of that mode. 
Because the two controls involved are not time-reverses 
of each other, the inequality of the two efficiencies is con- 
sistent with the time reversal arguments. As in the cavity 
case in paper I, the main reason for this deviation from 
the 7s = behavior is the dependence of the retrieval 
efficiency on the control. 

Finally, we discuss the effects of nonnegligible spin- 
wave decay on storage followed by retrieval considered 
in Sec. IVI CI Using the fact that nonzero 7 S keeps the 
general time reversal relationship between storage and 
retrieval exhibited in Eqs. (f3"5|) and (|3"6"|) . it is not hard to 
verify that nonzero 7 S still allows one to use time rever- 
sal iterations to optimize storage followed by retrieval. 
In particular, suppose that one is given a storage con- 
trol field and a (forward or backward) retrieval control 
field. Then one can find the optimal input mode to be 
used with these control fields by the following procedure: 
start with a trial input mode, store and retrieve it with 
the given pair of controls, time-reverse the whole proce- 
dure, and then repeat the full cycle until convergence is 
reached. Now suppose, on the other hand, one is given 
an input mode and is asked to choose the optimal storage 
and retrieval controls. Because of the spin wave decay, it 
is desirable to read out as fast as possible. As we discuss 
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in the next section, fast readout may be achieved in a 
time T ~ l/jd, so that if we assume that 7 S -C dry, the 
spin- wave decay during the retrieval will be negligible. 
If we further assume that the given input mode satis- 
fies the adiabatic limit Tdj 3> 1, then one should shape 
the storage control to store into the appropriate optimal 
spin-wave mode (Sd(l — 5) or a mode from Fig. O de- 
pending on the direction of retrieval) as if 7 S were zero 
and the input were proportional to £i n (t) exp(— 7 B (T — t)) 
(see Eq. (|44p). The total optimal efficiency will now de- 
pend on input pulse duration and shape: it will be mul- 
tiplied (and hence reduced relative to the 7 S = case) 

by dt\& m (t)\ 2 exp(-2%(f -t)). Finally, we note that 
when we consider storage followed by retrieval, in order to 
take into account the spin wave decay during the storage 
time [T, T r ], one should just multiply the total efficiency 
by exp(-2 7s (f r -f)). 

VII. FAST RETRIEVAL AND STORAGE 

We have shown that in the adiabatic limit (Td'y ^> 1, 
where T is the duration of the incoming pulse), one can 
optimally store a mode with any smooth shape and any 
detuning A. In this section, we solve Eq. (fT2|) analyt- 
ically in the second important limit, the "fast" limit, 
and demonstrate that this limit allows one to store opti- 
mally a certain class of input modes that have duration 
T ~ 1/ (d~f) . We also show that efficient (but not opti- 
mal) fast storage of any smooth pulse is possible as long 
as T7 < 1 and Tdj > 1. 

Exactly as in the cavity case in Sec. VI of paper I, in the 
fast limit, one assumes that f2 is very large during a short 
control pulse (|f2| 3> dj and |fi| 3> |A|) and keeps only 
terms containing (l on the right-hand side of Eqs. (J6]) and 
([7|) (or, equivalently, neglects all terms in Eq. (fT2|) except 
\(l\ 2 S and S). This gives Rabi oscillations between P and 
S and allows one to implement a fast storage scheme, in 
which the input pulse is resonant (A = 0) and the control 
pulse is a short 7r pulse at t = T, as well as fast retrieval, 
in which the control is a n pulse at t = T r . 

During fast retrieval, assuming that the 7r pulse is per- 
fect and that it enters the medium at t = (instead of 
t = T r ), the initial spin wave S = S(z) is mapped after 
the 7r pulse onto the optical polarization P — iS(z). We 
then solve Eq. (fTTj) for P(u,t), express £(u,t) in terms 
of P(u,t) using Eq. (fTT)|) . and take the inverse Laplace 
transform u — > z = 1 to arrive at 




When computing the fast retrieval efficiency, one can take 
the time integral analytically to find that the efficiency 
is again given by Eq. (fl"4"]) . which is consistent with the 
general proof in Sec. IIIII and the branching ratio argu- 
ment. In the cavity case discussion in paper I, we noted 
that the fast solution was a special case of the adiabatic 



solution with a suitable control. Similarly, the expression 
in Eq. (|4*6")l is also a special case of Eq. (f2"5| if we use 

n(t) = (l + iA)e~ lAi (47) 

and take the limit A — * 00 (although this violates the 
approximations made in deriving Eq. (|23|) ). 

Since the 7r-pulse control field in fast retrieval is fixed, 
optimal fast retrieval yields a single possible output 
mode, that of Eq. (|46[) with the optimal spin wave 
S(z) — Sd(z). By time reversal, the time-reversed version 
of this input mode (of duration T ~ 1/(7^)) is, therefore, 
the only mode that can be optimally stored using fast 
storage at this optical depth d. 

In order to confirm the time reversal argument and 
for the sake of completeness, one can also compute the 
optimal input mode for fast storage directly. For an input 
mode £i n (t) nonzero for t £ [0, T] and assuming that a 
perfect ir pulse arrives at t — T, we find that 

S(z,f) = -Vd dte-( f -V Jo (2 yjd(f-i)zj £ in (t) . 

(48) 

One can see that the fast retrieval and storage equations 
(|46|) and (|48|) obey, as expected, the same general time 
reversal relationship that we have already verified in the 
adiabatic limit in Eqs. ([55| and (j56"| . One can also ex- 
plicitly verify that the maximization of the storage ef- 
ficiency derived from Eq. (|48[) yields an optimal £i n (t) 
that is the normalized time reverse of Eq. (146|) evaluated 
with the optimal spin wave S(z) = Sd(z). It is worth 
noting that short exponentially varying pulses, reminis- 
cent of our optimal solution, have been proposed before 
to achieve efficient photon-echo based storage f34| . 

The solutions above give an incoming mode & m (t) that 
is optimal for fast storage alone or for fast storage fol- 
lowed by backward retrieval. Similarly, at each d, there 
is a mode that gives the optimal efficiency for fast storage 
followed by forward retrieval. This optimal input mode 
is the time-reverse of the output of fast forward retrieval 
from the spin-wave mode optimal for storage followed 
by forward retrieval (as computed through Eq. (|41[) and 
shown in Fig. [5]) . 

Finally, we note an important difference between fast 
storage in a cavity discussed in paper I and fast storage in 
free space. In a cavity, there is only one accessible spin- 
wave mode, and hence only one input mode that can be 
stored using fast storage (i.e., any input mode orthogo- 
nal to it will have zero storage efficiency). As shown in 
paper I, this input mode is exponentially rising with a 
time constant ~ 1/ (7C) , where C is the cooperativity 
parameter. Therefore, generating this mode, and hence 
obtaining high efficiency, may be hard in practice at high 
values of C. In contrast, in free space, any sufficiently 
smooth spin wave will have a high retrieval efficiency, 
and, by time-reversal, the time-reverses of the pulses fast 
retrieved from these spin waves can also be fast stored 
with high efficiency. One can, thus, explicitly verify using 
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Eq. (|48p. which allows one to compute these storage effi- 
ciencies, that if, in the original units, T7 <C 1 but at the 
same time Td'y ^> 1, the free-space fast-storage efficiency 
is close to unity. 



VIII. EFFECTS OF METASTABLE STATE 
NONDEGENERACY 

In the discussion of backward retrieval we have so far 
assumed that the two metastable states \g) and \s) are 
nearly degenerate. This has meant that during back- 
ward retrieval we could simply use the same equations 
as for forward retrieval but with the spin wave flipped: 
S(z) — ► S(L — z). If \g) and \s) are not degenerate and 
are split by to sg — cAfc, then during backward retrieval, 
instead of retrieving from S(L — z), we will have to re- 
define the slowly varying operators (see Eq. (|A10[) ) and 
retrieve from S(L — z) exp(— 2iAfcz), which significantly 
lowers the efficiency unless AkL <C yd. This condition 
on Afc can be understood based on the concept of the 
effective EIT window for the Fourier transform of the 
spin wave. As explained in Sec. IVI A 21 the width of this 
window is of order (in the original units) ~ yd/ L. The 
extra phase just shifts the Fourier transform off center 
by 2Afc, so that the efficiency will not be significantly 
affected provided the shift is much smaller than the win- 
dow width. We have confirmed numerically for S(z) = 1 
and for S(z) — \/3z that the AkL needed to decrease re- 
trieval efficiency by 50% from its AkL = value indeed 
scales as y/d (with proportionality constants sa 0.46 and 
« 0.67, respectively). 

There are two ways to understand physically why non- 
degeneracy of the metastable states ruins the backward 
retrieval efficiency. The first explanation, also noted in 
Ref. [35[ , comes from the fact that metastable state non- 
degeneracy breaks the momentum conservation on back- 
ward retrieval. During storage, momentum Ak is written 
onto the ensemble. Momentum conservation on back- 
ward retrieval, however, will require —Afc momentum in 
the spin wave. The second explanation comes from the 
fact that if Afc 7^ 0, then backward retrieval of optimal 
storage is no longer its time reverse. If we had not defined 
slowly varying operators, the spin wave that we store into 
our atoms would have had exp(jAfcz) phase written on 
it. Since time reversal consists of moving in the opposite 
direction and taking a complex conjugate, backward re- 
trieval will be the time-reverse of storage only if Afc = 0, 
in which case complex conjugation is trivial. Thus, if 
Afc 7^ 0, the optimization of storage does not simultane- 
ously optimize backward retrieval (unless, of course, we 
can apply the desired position-dependent phase to the 
atoms during the storage time, e.g., by a magnetic field 
gradient, or alternatively apply a n pulse that flips the 
two metastable states |35|V 

We would like now to optimize storage followed by 
backward retrieval in the presence of nondegeneracy 
(Afc 7^ 0). Following the general recipe of Sec. El in order 



to carry out the optimization, one has to start with an 
input pulse and a control pulse, do storage, then do back- 
ward retrieval with another control pulse. Then one has 
to time reverse the full process of storage and retrieval, 
and iterate till one gets convergence to a particular in- 
put (and spin wave). Specifically, we start with a trial 
spin wave Si (z) . To find the spin wave (in terms of op- 
erators that are slowly varying for forward propagation 
as defined in Eq. (|A10p ) that the optimal storage plus 
backward retrieval should use, we first rewrite Si(z) for 
backward-propagation slowly varying operators (i.e., add 
the 2Afcz phase), and then retrieve it backwards, time re- 
verse, and store. Using Eq. |2T|) . the iteration we get is 
(dropping an unimportant constant phase and going to 
our rescaled units) 

Sa(2)= / dz'fc r (z,5')e~ l2A ^' ' S{{z'), (49) 

where Afc = LAk. This iteration finds the eigenvector 
with the largest eigenvalue for the eigenvalue problem 

XS(z) = [ k x (z,z')e- i2A%B 'S*(z'). (50) 
Jo 

|A| 2 will then give the total maximum efficiency of storage 
followed by backward retrieval. In contrast to the Afc = 
case, the efficiencies of storage and retrieval in the opti- 
mal process are not generally equal. It is important to 
note that since the process we are optimizing followed by 
its time reverse corresponds to two iterations of Eq. (j49[) . 
after a sufficient number of steps, A settles into an os- 
cillation between |A| exp(ia) and |A|exp(— id) for some 
phase a. The eigenvector will oscillate between two val- 
ues differing only by an unimportant constant phase, so 
that either one can be used. 

While this procedure allows us to find the optimal spin 
waves, we should, for completeness, also determine, as in 
the Afc = case, which input fields the optimum can 
be achieved for. To do this, we, as before, consider the 
exactly solvable adiabatic and fast limits. In the adia- 
batic limit, the argument that retrieval can be shaped 
into any mode did not require the spin wave S(z) to be 
real, and it is therefore still applicable. By time reversal, 
we can, therefore, still achieve the maximum efficiency 
of storage followed by backward retrieval for any incom- 
ing mode of duration T such that Tdj 3> 1. Similarly, 
in the fast limit, using fast retrieval and time reversal 
we can find at each d a pulse shape with Td'j ~ 1 that 
gives the maximum efficiency. For completeness, we note 
that one can also generalize to Afc 7^ the method that 
uses the decayless mode s(z) to shape the optimal storage 
control, as described in Sec. IVI Bl However, since the op- 
timal control is unique, this method will, of course, yield 
the same control as the method based on retrieval and 
time reversal (as we showed explicitly in Sec. IVI Bl and 
Appendix [El for Afc = 0). We will, thus, omit here the 
extension of this method to the Afc 7^ case. 



19 



AkI 

20 r 



15 



10- 




100 



FIG. 7: If the two metastable states are not degenerate, the ef- 
ficiency of storage followed by backward retrieval will be low- 
ered relative to the degenerate case, because the energy dif- 
ference hu>sg introduces a momentum difference Afc = uj sg /c 
between the quantum and classical fields. As a function of 
d, the figure shows the momentum AkL, at which the opti- 
mal total efficiency of storage followed by backward retrieval 
falls to half of the AkL — value (dashed), and at which it 
is decreased to the optimal efficiency with forward retrieval 
(solid). 



We have, thus, demonstrated that we can optimize 
storage followed by backward retrieval for any given d and 
Afc. We also recall that we have shown in Sec. IVI Cl that 
for Afc = optimal storage followed by retrieval is accom- 
plished with backward retrieval. However, as we increase 
AfcL, the optimal total efficiency with backward retrieval 
will drop down to the optimal total efficiency with for- 
ward retrieval at some value of (AfcL)i. Increasing AfcL 
further up to another value (AfcL)2 will decrease the op- 
timal total efficiency with backward retrieval to half of its 
AkL = value (and then further down to zero). Figure 
0shows a plot of (AfcL)i (solid) and (AfcL)2 (dashed) as 
a function of d. As noted above, without reoptimization, 
(Afc_L)2 would go as y/d, but with optimization we see 
that it is linear in d, i.e., optimization makes the error 
less severe. (AfcL)i grows even slower than \fd. This 
is not surprising because at AfcL = optimal forward 
and optimal backward errors both fall off as 1/d, except 
with different coefficients and, thus, eventually get very 
close to each other, so it takes a small AkL to make them 
equal. 

In Figs. [5Ja) and[5{b), we show the magnitude 1-5(5)1 
and the phase Arg[S'(z)], respectively, of the optimal 
mode (defined for the forward-propagating slowly vary- 
ing operators as in Eq. (|A10|) ) at d — 20 for different 
values of AfcL. As we increase AfcL, the optimal mode 
becomes concentrated more and more near the back end, 
i.e., it becomes favorable to effectively decrease the opti- 
cal depth (i.e., decrease effective L) in order to decrease 
effective AfcL. The phase of the optimal mode is approx- 



imately linear, i.e., S(z) oc exp(— ik^z) for some fco- At 
AfcL = 0, fco = 0. Interestingly, instead of just grow- 
ing from linearly with AfcL, fco first increases but then 
above AfcL ~ 7.5 starts decreasing again. 

|S(2)| 




0.25 

Arg[S(Z)] 




FIG. 8: (a) The magnitude and (b) the phase of the optimal 
mode for storage followed by backward retrieval at d — 20 
for the indicated values of AkL. The phase of the optimal 
mode at AkL = is 0. The phase is plotted for the forward- 
propagation slowly varying operators as defined in Eq. (|A10|l . 



IX. SUMMARY 

In conclusion, in this paper, we have presented a de- 
tailed analysis of the storage and retrieval of photons in 
homogeneously broadened A-type atomic media in free 
space and made a comparison to the cavity model de- 
scribed in paper I. From the investigation in the present 
paper emerges a new physical picture of the process of 
storage and retrieval in this system: first of all, the 
retrieval is essentially an interference effect where the 
emission from all the atoms interferes constructively in 
the forward direction. This constructive interference en- 
hances the effective decay rate into the forward direc- 
tion to dj. The branching ratio between the desired for- 
ward radiation and the unwanted spontaneous emission 
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is then simply given by the ratio between the various de- 
cay rates and is 77 ~ ^7/(7 + dry) ~ 1 — l/rf, irrespective 
of the method being used to drive the excitation out of 
the atoms. Secondly, the storage process is most conve- 
niently viewed as the time reverse of retrieval. 

In the present paper, we have used this physical picture 
to derive the optimal strategy for storage and retrieval 
and the optimal efficiency that is independent of whether 
one works in the Raman, EIT, photon-echo, or any other 
intermediate regime. In particular, we showed how to 
achieve the optimal storage of any smooth input mode 
at any detuning of duration T 3> 1/(07) (the adiabatic 
limit, including Raman and EIT) and of a particular class 
of resonant input modes of duration T ~ 1/(07) (t ne 
fast or photon-echo limit). This analysis is extendable to 
other systems. In particular, in paper III, we consider the 
effects of inhomogeneous broadening on photon storage 
in A-type atomic media. Extensions to other systems, 
such as the double-A system [36| or the tripod system 
[3^ . should also be possible. 

We also suggested a novel time reversal based itera- 
tive procedure for optimizing quantum state mappings. 
Moreover, we showed that for the case of photon storage, 
this procedure is not only a convenient mathematical tool 
but is also a readily accessible experimental technique 
for finding optimal spin waves and optimal input-control 
pairs: one just has to be able to measure the output mode 
and to generate its time reverse. Following the present 
work, this procedure has recently been implemented ex- 
perimentally with classical light We also expect this 
optimization procedure to be applicable to other systems 
used for light storage, such as tunable photonic crystals 

The presented optimization of the storage and retrieval 
processes leads to a substantial increase in the memory 
efficiency whenever reasonable synchronization between 
the input photon wave packet and the control pulse can 
be achieved. We, therefore, expect this work to be impor- 
tant in improving the efficiencies in current experiments, 
where optical depth is limited by various experimental 
imperfections such as a limited number of atoms in a trap 
[39( , competing four- wave mixing processes in a warm va- 
por cell f^L or inhomogeneous broadening in solid state 
samples |40| . 



X. ACKNOWLEDGMENTS 



We thank M. Fleischhauer, M. D. Eisaman, E. Polzik, 
J. H. Muller, A. Peng, J. Nunn, I. Novikova, 
D. F. Phillips, R. L. Walsworth, M. Hohensee, M. Klein, 
Y. Xiao, N. Khaneja, A. S. Zibrov, P. Walther, and 
A. Nemiroski for fruitful discussions. This work was sup- 
ported by the NSF, Danish Natural Science Research 
Council, DARPA, Harvard-MIT CUA, and Sloan and 
Packard Foundations. 



APPENDIX A: DETAILS OF THE MODEL AND 
DERIVATION OF THE EQUATIONS OF 
MOTION 

In Sec. HH we presented a short introduction to the 
model and stated the equations of motion without deriva- 
tion. In this Appendix, we provide the details of the 
model, as well as the derivation of the equations of mo- 
tion ([5])-([7]). Since the model and the assumptions made 
are very similar to those presented in the cavity case in 
paper I, we will often review some of them only briefly. 

The electric field vector operator for the quantum field 
is given by [4l[ 



Ei(z) = £i 



\ 47rce< 



\l/2 , 

^- / dua u e iu */ c + h.c, (Al) 
qA/ J 



where h.c. stands for Hermitian conjugate and where we 
have a continuum of annihilation operators a w for the 
field modes of different frequencies lu that satisfy the com- 
mutation relation 



(A2) 



By assumption, the field modes corresponding to a u 
for different lu have the same transverse profile and are 
nonempty only around u> = u>\. We have here assumed 
that the cross section A of the beam is identical to the 
cross section of the ensemble. In typical experiments, the 
beam is smaller than the size of the ensemble, and in this 
case the relevant number of atoms N should only be the 
number of atoms interacting with the beam. However, as 
we see from the final equations ©-([T]), the only relevant 
quantity is the optical depth d, which does not depend 
on the area A, so that when everything is expressed in 
terms of d, the precise definition of N and A is irrelevant 
(see the end of this Appendix). 

The copropagating classical control field vector 



E2OM) = £ 2^2(< - z/c) cos(w 2 (i - z/c)) 



(A3) 



is a plane wave with polarization unit vector e 2 and car- 
rier frequency a>2 modulated by an envelope £ 2 {t — z/c), 
which we assume to be propagating with group velocity 
equal to the speed of light c since almost all the atoms 
are assumed to be in the ground state \g) and are, thus, 
unable to significantly alter the propagation of a strong 
classical field coupled to the |s)-|e) transition. 

The Hamiltonian in Eqs. (A3)-(A5) in paper I is then 
modified to 



H = Ho + V, 
H 



(A4) 



r N 

i—1 



N 



-iLU 2 (t-Zi/c) 



(A6) 
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Here a* = \^)a(v\ is the internal state operator of the 
ith atom between states \fi) and is the position of 

the ith atom, d is the dipole moment vector operator, 
tt{t-z/c) = (e\(d-e 2 )\s)£ 2 (t- z/c)/(2K) is the Rabi fre- 
quency of the classical field, and g = (e|(d-ei)|<?) v/i^f^CE 

(assumed to be real for simplicity) is the coupling con- 
stant between the atoms and the quantized field mode, 
where we have chosen the length of the quantization vol- 
ume to be identical to the ensemble length (this choice 
does not affect the results obtained below). We note 
that in order to avoid carrying extra factors of 2 around, 
£1 is defined as half of the traditional definition of the 
Rabi frequency, so that a it pulse, for example, takes 
time tt/(2Q). 

Since the position dependence along the ensemble mat- 
ters, we divide our ensemble into thin slices along the 
length L of the ensemble (z = to z = L) and introduce 
slowly varying operators 



z i=l 
1 S 

i S 



iu>2(t — Zi/c) 



■iuj-L [t — Zij c) 



N 



(A7) 
(A8) 
(A9) 



&sg(z,t) 

£(z,t) 



where sums are over all N z atoms in a slice of atoms 
positioned at z that is thick enough to contain N z 1 
atoms but thin enough that the resulting collective fields 
can be considered continuous. The normalization of £ is 
chosen to ensure that it is dimensionless, which will be 
necessary to yield the simple dimensionless expressions 
in Eqs. I©-®. For these slowly varying operators, the 
effective Hamiltonian is 



i^£^) e ~ i(wl ~ W2)(t " </c) > ( A1 °) 



iu, x (t-z/c) / dua^ty^^Ml) 




H = 



J dujhjja^a^ — huji— J dz£^ (z,t)£(z,t) 



+ / dzhn(z) 



A 



a ee (z,t) - (Cl(t- z/c)a es (z,t) 



+g£(z,t)a eg (z,t) + h.c. 



' eg 

and the same-time commutation relations are 

1 



(A12) 



[a^{z,t),a al 3{z\t)] 



n(z) 

xd(z-z') 
£(z,t),&(z',t)] =LS(z-z'). 



M)) 

(A13) 
(A14) 



Under the same assumptions as in the cavity case in pa- 
per I and defining P = \f~Na ge and S — ^/N& gs , the 
Heisenberg equations of motion yield Eqs. H])-©, where, 
as in the cavity case in paper I, 7 may include extra de- 
phasing in addition to radiative decay. We note that P 
and S are defined to be dimensionless in order to yield 
fully dimensionless Eqs. ©-{7]). The V77 in the defi- 
nitions of P and S is required in order to have the fi- 
nal dimensionless equations depend on g, N, and L only 
through the optical depth d. Similarly to the cavity case 
in paper I, from the generalized Einstein relations, the 
only nonzero noise correlations are [42| 



(F P (z,t)Fl(z>,t>)) = 



N 
n(z) 

N 
n(z) 



5(z - z')6(t - t'),(A15) 



(F s (z,t)Fl(z',t')) = -^S(z -z')5(t-t'). (A16) 



Again the fact that normally ordered correlations are 
zero, as in the cavity case in paper I, means that the 
incoming noise is vacuum, which is precisely the reason 
why, as noted in Sec. II in paper I, efficiency is the only 
number we need in order to fully characterize the map- 
ping. The property of our system that guarantees that 
the incoming noise is vacuum is the absence of decay out 
of state \g) into states |e) and |s). We refer the reader to 
Appendix A of paper I for a detailed discussion of why 
this is a reasonable assumption in most experimental re- 
alizations. 

We will now show how our field and atomic operators 
can be expanded in terms of modes, which is necessary in 
order to obtain and interpret the final complex number 
equations ©-([7]). Under the assumption that almost all 
atoms are in the ground state at all times, commutation 
relations (|A13p imply 



S(z,t),S\z',t) 
P(z,t),PHz\t) 



N 
n(z) 
N 
n(z) 



S(z-z'), 
S(z-z'). 



(A17) 
(A18) 



Equation (|A17|) allows us to expand S(z, t) in terms of 
any basis set of spatial modes {g a (z)} satisfying the or- 
thonormality relation dzg^(z)gfj(z) = S a fj and the 
completeness relation ^ Q g* a (z)g a (z') = 8{z — z') as 



S{z,t) 



I N 
n(z) 



Y,g a (z)c a {t), (A19) 



where the annihilation operators {c Q } for the spin- wave 
modes satisfy 



c a (t) , dp(t) =5 a 0- 



(A20) 



For the freely propagating input field £ m {t) = £(0,t) 
and output field £ ut(t) = £(L,t) we have the following 
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commutation relations: 



£„(*), 4(0 



L 

c 
L 



S(t-t'), 



(A21) 



£out(*),£ T ut(*') = (A22) 
J c 

which differ from their cavity case counterparts in 
Eq. (A19) of paper I only in normalization. These com- 
mutation relations allow us to expand, as in the cavity 
case in paper I, the input and the output field in terms 
of any basis set of field (envelope) modes {h a (t)} de- 
fined for t 6 [0,oo), satisfying the orthonormality rela- 
tion J dth* a (t)hf3(t) = 5 a f3 and the completeness rela- 
tion E Q K (t)h a (f) = S(t - t'), as 

4x(t) =\h;Yl M*) a «> ( A23 ) 



(A24) 



^out(i) = \ — y^h a (t)b ai 



where annihilation operators {a a } and |^q| for the in- 
put and the output photon modes, respectively, satisfy 
the usual bosonic commutation relations (see Eq. (A22) 
in paper I). 

All atoms are initially pumped into the ground state, 
i.e., no P or S excitations are present in the atoms. We 
also assume that the only input field excitations initially 
present are in the quantum field mode with annihilation 
operator do corresponding to an envelope shape ho(t) 
nonzero on [0, T]. Precisely as in the cavity case in paper 
I, the only parts of the operators that will contribute to 
the efficiency will be the parts proportional to do. We can 
therefore reduce our problem to complex number equa- 
tions. These equations and the corresponding initial and 
boundary conditions are given in Sec.|TTl To get back the 
nonvacuum part of the original operator from its complex 
number counterpart, one can just multiply the complex 
number version by do. 

We conclude this Appendix with a verification that 
d = g 2 NL/(-yc) is independent of the size of the beam 
and, for a given transition, depends only on the density 
of the atoms and the length of the ensemble. This can 
be seen directly by inserting the definition of g into d 
and defining the atomic number density p(z) = n{z)/A 
(which, by assumption, is uniform in the direction trans- 
verse to the propagation direction). The expression for 
d then becomes d = |(e|(d • £i) |<7) | 2 <^i J p{z)dz / (2/ieo7c). 
So it is indeed independent of the size of the beam and, 
for a given transition, only depends on the density and 
length of the ensemble. 



atomic excitations (i.e., the retrieval is complete), the re- 
trieval efficiency rj r depends only on the optical depth d 
and the spin wave S{z) but not on the detuning A and 
the control field envelope SI (t). In this Appendix, we 
show that for complete retrieval, not only the total effi- 
ciency but also the distribution of spontaneous emission 
loss (or more precisely loss due to polarization decay 7) 
as a function of position is independent of the control 
and the detuning. 

Equations of motion (O-® imply that 

d- z \£{i,i)\ 2 + d i \P(^t)\ 2 + d i \S(S,t)\ 2 = -2\P(S,t)\ 2 . 

(Bl) 

Integrating both sides with respect to z from to 1 and 
with respect to t from T r to 00, using the initial conditions 
S(z,f r ) = S(z) (where /* dz \S(z) | 2 = 1) and P{S,f r ) = 
0, the boundary condition £ (0, t) = 0, and the complete 
retrieval condition S(z, 00) = P(z, 00) = 0, we find, using 
Eq. ©, that 



dz l(z), 



(B2) 



where the position-dependent loss per unit length is 



Computing l(z), we find 



2 / dt\P{z,t)\ 2 . 




(B3) 



l(z) = 2£- 



dtP(u,t) [P{u'*,t)]* 



-S(u) [S(u'*)Y 



,(B4) 



where C^ 1 with subscript u, u' — > z means that inverse 
Laplace transforms are taken with respect to u and u' and 
are both evaluated at z. In the last equality, Eq. (11311 and 
the conditions at t — T r and t — 00 were used. Therefore, 
we see that l(z) is independent of the detuning and the 
control. Moreover, the inverse Laplace transforms 
can be taken analytically to give 



-■i|2 



l(~z) = \S(z)\ 



Re 



S(z)d / dzS*(z — z)e 




dz' j dz"S{Z-2')S*{2-2")—e 

2io (VFf) - z -±f h (dVW) 



:«\ d2 „-f (f'+z") 



(B5) 



APPENDIX B: POSITION DEPENDENCE OF 
LOSS 



APPENDIX C: IMPLEMENTATION OF THE 
INVERSE PROPAGATOR USING TIME 
REVERSAL 



We have shown in Sec. IIIII that, provided the retrieval 
control pulse is long and/or powerful enough to leave no 



In Sec. HVi we exploited the fact that time re- 
versal could be used to realize the inverse evolution 
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U^ 1 [T, 0; fl(t)]. In this Appendix, we would like to ex- 
plain carefully what we mean by the time reversal oper- 
ator T and to prove Eq. (|19p . 

To define the time reversal operator T, we first choose 
a basis for our single-excitation Hilbert space consist- 
ing of {<r* g (ground)}, { a\ g \ ground)}, and {a\ |ground)}, 
where | ground) is the state with no photons and no 
atomic excitations (i.e., all atoms in the ground state), 
i runs over all atoms, z runs over all positions, and 
a\ = (27rc) -1 / 2 J duj exp (— iuz/c) oj,. We then define 
the time reversal operator T (equivalent to the complex 
conjugation operator K in Ref. [43j) as follows: T\ip(t)) 
means taking the complex conjugates of the expansion 
coefficients of a state \i/j(t)) in the above basis, while 
TOT means taking complex conjugates of the matrix 
elements of the operator O when O is written in the 
above basis (we will, thus, write TOT = O*). For ex- 
ample, this definition implies that in addition to com- 
plex conjugating the envelope of the photon, time re- 
versal flips the photon momentum: TajT = a_ w and 
Ta u |ground) = d_o,|ground). Some of the properties of 

f are f 2 = 1 and |(Vi|t|V 2 )| = \{f^ 2 )\- 

We now turn to the proof of Eq. (Tl9|) . We start by 
noting that 

U- 1 [T,0;Q{t)\ = U[0,T;Cl(t)]=ffU[0,T;Cl(t)]ff 
= TU*[0,T;n(t)]f (CI) 

and, therefore, using Eq. (jT5J) , 

V = \{b\U[T,0;£l(t)]\a)\ 2 = \(a\fU*[0,T;Q.(t)}f\b)\ 2 . 

(C2) 

To evaluate U*[0,T;fl(t)] and to find a way to imple- 
ment it physically, let us first consider the simplest case, 
where the Hamiltonian responsible for the evolution is 
independent of time and respects time reversal symme- 
try THT = H. This is equivalent to H* = H. In 
this case, the evolution operator is given by U[T, 0] = 
exp(-iHT/h), and, therefore, U*[0,T] = U[T,0}. So if 
the Hamiltonian obeys time reversal symmetry, one can 
physically implement U* [0, T] simply by evolving the sys- 
tem for a time T. Applied to Eq. (|C2[) . this would mean 
that the probability to go from \a) to \b) due to unitary 
evolution U is the same as the probability to make the 
transition from T\b) to T\a) due to this evolution. In 
other words, according to time reversal, if our Hamilto- 
nian obeyed time reversal symmetry, we would be able 
to map the time reverse £*(T — t) of the output mode 
onto the spin wave S* with the overlap storage efficiency 
equal to the retrieval efficiency from the spin wave S. 

In general, the interaction does not obey time reversal 
symmetry because of the classical control, which may 
depend on time and may be complex. To extend the 
discussion to this situation, we shall use the equation of 
motion for the propagator 

. n dU[r u rr,m] = H[n;mP[ri,r 2 ;Sl(t)l (C3) 



where we have highlighted the dependence of the Hamil- 
tonian on the control field by including Q(t) as an argu- 
ment of H . By taking the complex conjugate of Eq. (|C3|) . 
we can turn it into the equation for the time-reversed in- 
verse propagator ftj- 1 [r 2 , n; Q(t)]f = U*[n, r 2 ;ft(t)]: 

Note that we have not made any assumptions about time 
reversal symmetry being a symmetry for the system. 

In order to implement the evolution in the last ex- 
pression of Eq. (|C2[) . we need to realize a time evolution 
U[T, 0; fi'(t)] such that U[T, 0; fi'(t)] = U*[0,T,Q(t)]. To 
do this, we consider operators U[t, 0; ft'(t)] and U*[T — 
r, T, Sl(i)] whose equality at time r = T will imply the 
desired evolution. The equation of motion for these op- 
erators are, using Eqs. (|C3|) and (1C4[) . 

. n dU*[T-T/T;a(t)} = H*[T- T ;n{t)}U*[T-T,T;Sl(t)], 

. h dU[T,0;n'(t)] = £ n'{t)]U[r, 0; Cl'(t)]. (C5) 
ar 

Since the two operators are the same at time r = 0, 
U[0, 0, Q'(t)} = L>*[T, T, n(t)] = i, the two operators will 
be identical (and, in particular, equal at r = T) if they 
obey the same differential equation, which is the case if 

H[t,n'(t)]=H*[T-t;fl{t)]- ( C6 ) 

We would like to evaluate the right-hand side of Eq. (|C6|) . 
By inspecting the Hamiltonian in Eqs. (|A4|) - (|A6[) . we 
see that the only non-trivial parts are the classical con- 
trol field (including the carrier LO2) and the exponential 
in the interaction with the quantum field (the part of 
Eq. (|A6|) containing g) . According to the right hand side 
of Eq. (|C6p . we would like to apply the complex conju- 
gation to H and evaluate it at time T — t. The interac- 
tion with the quantum field is actually invariant not only 
under the change of time but also under complex con- 
jugation since applying time reversal to a u e mz, ' c simply 
changes ui to —u>, which can be changed back to ui by flip- 
ping the sign of the integration variable. The application 
of complex conjugation and the change of time to the 
interaction with the classical field is equivalent to using 
the time- reversed control field envelope Sl'(<) = f2*(T — t) 
and a carrier wave vector propagating in the opposite di- 
rection. Combining this result with Eq. (|Cip . we arrive 
at Eq. (fill)) . This means that by using the time-reversed 
control field f2*(T — t), we can map T\b) onto T\a) with 
the probability equal to the probability that \a) goes to 
1 6) using £l(t). 
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APPENDIX D: PROOF OF CONVERGENCE OF 
OPTIMIZATION ITERATIONS TO THE 
OPTIMUM 

In Sec. [V] we omitted the proof that iterative applica- 
tion of NPaU -1 PbU to a unit vector \a) G A converges 
to |a max ) (unless (a|a max ) = 0) and that NPBU\a ma , x ) 
optimizes U^ 1 as a map from B to A. In this appendix, 
we present this proof. 

For any unit vector |a) in the subspace A of "ini- 
tial" states, the efficiency is defined as rj = \PbU\o)\ 2 , 
where Pb denotes the projection on the subspace B 
of "final" states. We are looking for |a max ) € A that 
gives the maximum efficiency ry max = \PBU\amax)\ 2 ■ 
By appropriately adjusting the phases, we can write 
y/Vmax = (&max|^|a max ) for some unit vector |6 max ) = 
AfPBU\a max ) E B. From the unitarity of U (see, 
for example, Eq. (flBf ). it follows that |& m ax) opti- 
mizes as a mapping from B to A. We will 
prove two claims, from which the desired convergence 
result will follow immediately. Claim (1): if \a\) is 
orthogonal to |a max ), U\a\) is orthogonal to |6 max ). 
Proof: suppose (6 ma x|^|ai) = ^ 0, then defining 
\a) = (^?7max|a max ) + /3*|ai))/\/?7max + \P\ 2 , we have 
(b max \U\a) = y/r) max + \P\ 2 > y/r] max , which contradicts 
the fact that |a max ) was optimal. A similar proof can 
be given for claim (2): if \b%) is orthogonal to |6 max ), 
is orthogonal to |a max ). From these two claims it 
immediately follows that if we start with a state orthog- 
onal to |a max ), we will never approach |a max ). On the 
other hand, we will show now that if we start with \a) = 

^l^max ) + V 1 - M 2 K) ( for some unit K> orthogonal 
to |o max ) and for some a ^ 0), we will indeed approach 

). We have P B U\a) = a(b max \U\ 
y/l — |a| 2 (6i|C/|ai)|6i) for some unit \bi) GB. By claim 
(1), the two parts of P B U\a,) are orthogonal and, since 
Kiax) is optimal, |{& max |^|a max )| > |(fei|t/|oi)|. Thus, 
the fraction of |6 max ) in PbU\cl) is greater than the frac- 
tion of |a max ) in |o). After the application of PaU -1 and 
during subsequent iterations, the optimal fraction will 
similarly grow. This shows that we will indeed reach the 
optimum, unless we start with something orthogonal to 
it. 



APPENDIX E: SHAPING THE CONTROL FIELD 
FOR THE OPTIMAL ADIABATIC STORAGE 

In this Appendix, we use Eq. (pT7|) to find the control 
field for the storage of any given mode £- m (t) into any 
given decayless spin- wave mode s(z). We then verify that 
the optimal storage control found through this procedure 
using the optimal decayless spin-wave mode s(z) gives 
storage into 5^(1 — z), the optimal mode for backward 
retrieval, with efficiency equal to the optimal retrieval 



efficiency r/™ ax . We also verify that this control field is the 
time-reverse of the control field that retrieves the optimal 
spin- wave mode into £* n (T — t), that is, the time-reverse 
of the input mode. 

In order to solve for f2(t) from Eq. (|37[) . we note that 
q(z, t) satisfies 



dzq[z,i)q*(z,i') = 6(t-i'), 



dtq(z, t)q* (z , t) = S(z — z), 



(El) 



(E2) 



where we have used the identity 



1 



dxJo(ax)Jo(bx)x = —S(a — b), (E3) 



and where Eq. (|E2|) requires h(0,T) = oo (we will dis- 
cuss below that this requirement can be relaxed without 
significant loss in efficiency). Using Eq. (|E1[) . we see that 
J Q dz\s(z)\ 2 = 1, as expected from unitarity (since we 
neglect both the leakage and the decay rate 7 in Eq. (j37]) . 
the transformation between £i n (t) and s(z) is unitary). 
Using Eq. (|E1[) . we can also invert Eq. (|3T[) to get 



fin® 



dzq*{z, t)s(z) 



(E4) 



Clearly, Eqs. (|37|) and (|E4|) establish a 1-to-l correspon- 
dence, for a given Q(t), between input modes £ ln (t) and 
decayless modes s(z). This 1-to-l correspondence is the 
demonstration of the unitarity of the map defined by 
Eq. (I37p . For the purposes of shaping the control field, 
it is crucial that Eqs. (|37|) and (|E4j) also establish a 1- 
to-1 correspondence, for a given £; n (i), between controls 
ri(t) (satisfying h(0,T) = 00 and nonzero whenever f in (t) 
is nonzero) and normalized decayless propagation modes 
s(z). In particular, Eq. ([3T)l itself allows to determine 
s(z) from Q(t) and £- m (t). To solve for Cl(t) given s(z) 
and £i n (t), we integrate, as in Sec. IVI A 31 from to t 
the norm squared of both sides of Eq. (IE4|) and change 
the integration variable from t' to h' — h(t',T) on the 
right-hand side to obtain 



dt'\£ in (t' 



h(0,T) , 

dh's- 

A 2 



h(t,T) 

dze 



■■J„\hl!§)s(z) 



(E5) 



Using h(T, T) = and the normalization of s(z) and 
£in(t), the evaluation of Eq. (|E5j) at t = T, implies that 
h(0,T) — 00 (unless the expression inside the absolute 
value sign vanishes for all b! greater than some finite 
value). This divergence is, however, just a mathematical 
convenience: truncating |f2(f)| does not significantly af- 
fect the efficiency, as we discuss in Sees. IVI A 3[ IVI Bl and 
IVIDI We also give in Sec. IVI Bl approximate expressions 
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for how big the optimal |f2| is in the Raman and reso- 
nant limits. Replacing h(0,T) with oo, we use Eq. (|E5|) 
to solve numerically for h(t, T), exactly as in Sec. IVI A 31 
| Cl(t) | is then deduced by taking the square root of the 
negative of the derivative of h(t,T). The phase of Cl is 
found by inserting |Q| into Eq. (|E4|) is given by 



Arg 



+Arg 



f+Arg M ]+^ 



.(E6) 



The various terms in Eq. (|E6[) can be interpreted in a way 
similar to the terms in the phase of the retrieval control 
in Eq. ([3"3"[) in Sec. IVI A 31 The only minor difference in 
the interpretation is that in the limit A — > 0, the third 
term seems to diverge. However, one can check that in 
this limit the last term cancels (up to a constant) with 
the third term to ensure that the phase of the optimal 
control is still given solely by the phase of the desired 
output, as expected for the resonant limit. Finally, we 
note that the same remarks as at the end of Sec. IVI A 31 
regarding the ability to truncate the divergences of |0(t)| 
without significant loss in efficiency apply. 

Using Eq. (I39|) . we show in Sec. IVI B I how to find the 
optimal decayless mode s(z), which should then be used 
in Eq. ([37]) to shape the optimal control (using Eqs. (|E5|) 
and (|E6p ). Having derived the optimal storage control in 
this way, we can now explicitly verify the results obtained 
from the time reversal reasoning. Two of these results are 
that the mode S(z, T) used in optimal storage is just the 
optimal mode for backward retrieval and that the opti- 
mal storage efficiency and the optimal retrieval efficiency 
are equal. One can explicitly verify these statements by 
checking that the application of Eq. (|39f to the iteration 
used to find the optimal s(z) gives the iteration used to 
find the optimal backward retrieval mode, i.e., 



dz k T (z, z')Si(z'), 



(E7) 



where k T (z,z') is defined in Eq. [[IS]). Thus, S(z,f ) that 
we compute from the optimal s(z) using Eq. ([3"9"]> is in- 
deed the optimal mode for backward retrieval, and the 
optimal storage efficiency is indeed equal to the optimal 
retrieval efficiency r/^ 1 ^: 



S(z,T) = ^j™*S d (l-z). 



(E8) 



Another consequence of time reversal is the fact that 
the optimal storage control for a given input mode is 
the time reverse of the control that gives optimal back- 
ward retrieval into the time reverse of that input mode. 
One can verify this directly by comparing the expressions 
of the magnitude and phase of the two controls. How- 
ever, a simpler approach is to consider £ on t(t) given by 
Eq. ([23]) for the case of retrieval with a certain control 
fl(t) from the optimal mode, i.e., S(z) — S d (z). We then 
use Eq. (pM]) to compute the spin wave S2 (z) that results 
from storing £* ut (f - t) with Q*(f - t). In the limit 
T — > 00 (to make sure that we fully retrieve £ ou t before 
sending its time reverse back in), we can take the time 
integral explicitly to find 



S 3 (z) 



dz'k T (z,z')S d (l-z')=r)™*S d (l-z), (E9) 



where in the last step we used the definition of S d (l — z) 
as the eigenvector of fc r with the largest eigenvalue ry™ ax 
equal to the optimal retrieval efficiency. Thus, the total 
efficiency of optimal retrieval followed by time-reversed 
storage is (i~if lax ) 2 . So we have shown explicitly that the 
time reverse of optimal retrieval gives storage into S d 
with the maximum efficiency ?y" lax , confirming what we 
have shown in Sec. IIVI based on general time reversal ar- 
guments. Since we have shown in this Appendix that the 
optimal control field is unique, we have therefore con- 
firmed that the control that optimally stores a given in- 
put and the control that optimal retrieves into the time 
reverse of that input are time-reverses of each other. 
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